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Ideal (mono-domain) smectic- A elastomers crosslinked in the smectic- A phase are simply uniaxial 
rubbers, provided deformations are small. From these materials smectic-C elastomers are produced 
by a cooling through the smectic-y4 to smectic-C phase transition. At least in principle, biaxial 
smectic elastomers could also be produced via cooling from the smectic-A to a biaxial smectic 
phase. These phase transitions, respectively from Dooh to C2/1 and from -Dooh to Dih symmetry, 
spontaneously break the rotational symmetry in the smectic planes. We study the above transitions 
and the elasticity of the smectic-C and biaxial phases in three different but related models: Landau- 
like phenomenological models as functions of the Cauchy-Saint-Laurent strain tensor for both the 
biaxial and the smectic-C phases and a detailed model, including contributions from the elastic 
network, smectic layer compression, and smectic-C tilt for the smectic-C phase as a function of both 
strain and the c-director. We show that the emergent phases exhibit soft elasticity characterized by 
the vanishing of certain elastic moduli. We analyze in some detail the role of spontaneous symmetry 
breaking as the origin of soft elasticity and we discuss different manifestations of softness like the 
absence of restoring forces under certain shears and extensional strains. 

PACS numbers: 83.80.Va, 61.30.-v, 42.70.Df 



I. INTRODUCTION 

Liquid crystalline elastomers [lj are fascinating hybrid 
materials that combine the elastic properties of rubber 
with the orientational and positional order of liquid crys- 
tals 0- As in conventional liquid crystals, there exists 
a great variety of phases in liquid crystalline elastomers. 
For example, nematic, cholesteric, smectic- A (SnxA), chi- 
ral smectic- A* (SmA*), smectic-C (SmC), and chiral 
smectic-C* (SmC*) phases have been created in elas- 
tomeric forms [lj]. Among these elastomers, nematics 
have to date received the most attention leading to the 
discovery of a number of remarkable properties of these 
materials such as soft elasticity [1, 0, H,H, 0] > dynamic soft 
elasticity [1, H, El] , and anomalous elasticity [III El, EH ■ 
In contrast, the understanding of smectic elastomers is 
much less developed, at least from a theoretical point of 
view. Given that there is a substantial literature treat- 
ing the synthesis and experimental properties of smectic 
elastomers 0, El, El, E3, El, El, HI, El, HI , that smec- 
tic elastomers have intriguing properties and potential 
for device applications such as manometer-scale actua- 
tors [11| , and that smectics play a leading role in conven- 
tional liquid crystals where they have attracted outstand- 
ing scientific and technological interest since the discov- 
ery of spontaneous ferroelectricity in C* smectics [24j . it 
is clear that a deeper theoretical understanding of smec- 
tic elastomers is desirable. 

To date there have been, as far as we know, only rela- 
tively few theoretical investigations of smectic elastomers 
for the apparent reason that such investigations are dif- 
ficult due to the complexity and low symmetry of the 
material. Terentjev and Warner developed expressions 



for the elastic energies of SmA and SmA* elastomers [11] 
as well as for SmC and SmC* elastomers [26( based on 
group theoretical arguments. The coupling of the smec- 
tic layers to the elastic network was critically discussed 
and shortcomings of Ref. [Hj] in this respect were cor- 
rected in (27[. Subsequently, the damping effect of the 
rubber-elastic matrix on the fluctuations of the smectic 
layers, leading to a suppression of the Landau-Peierls in- 
stability and true one-dimensional long-range order, was 
analyzed in (28|. Weilepp and Brand [11] discussed an 
undulation instability as a possible explanation for the 
turbidity of SmA elastomers under stretch along the nor- 
mal of the smectic layers. Osborne and Terentjev [ll] 
derived expressions for the effective elastic constants of 
SmA elastomers when these are viewed as effectively uni- 
axial systems and revisited the suppression of fluctua- 
tions of the smectic layers. Very recently, Adams and 
Warner (AW) [U, [H, [33| set up a model for the elas- 
ticity of smectic elastomers by extending the so-called 
neoclassic model of rubber elasticity, which was origi- 
nally developed and very successfully used to describe 
nematic elastomers [lj, to include the effects of smec- 
tic layering. Also just recently we worked out theories 
for the low- frequency long- wavelength dynamics of smec- 
tic elastomers in respectively the SmA, biaxial and SmC 
phases [34| . 

In principle, smectic elastomers can be produced ei- 
ther by cooling nematic elastomers through the transi- 
tion to the smectic phase or by crosslinking in the smectic 
phase. Both SmA and SmC elastomers have been pre- 
pared mostly by the second method, e.g., by crosslink- 
ing side chain liquid crystalline polymers, which have a 
tendency to form layers because their mesogens are of- 
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ten immiscible with the polymeric backbone [ljj], or by 
crosslinking polymer chains or hydrophobic tails in bi- 
layer lamellar phases of respectively diblock copolymers 
or surfactant molecules [l5| . Crosslinking in the smectic 
phase tends to lock the smectic layers to the crosslinkcd 
network [27| • Without this lock-in, the phase of the smec- 
tic mass-density-wave can translate freely relative to the 
elastomer as it can in smectics in aerogels [Hj]. To keep 
our discussion as simple as possible, we will not consider 
in the following the case of crosslinking in the nematic 
phase, and we take the lock-in of the smectic layers and 
the elastic matrix as given. 

As in nematics, unconventional properties are most 
pronounced in samples of smectic elastomers that have 
an ideal, monodomain morphology so that the director 
has a uniform orientation throughout. In practice, how- 
ever, liquid crystal elastomers tend to be non-ideal, i.e. 
the material segregates into many domains, each having 
its own local director. In order to avoid such polydomain 
samples, elaborate crosslinking schemes involving elec- 
tric or mechanical external fields for aligning the director 
have been developed and successfully applied to smectic 
elastomers 0, [TJ Gil- The latest achievement in this 
respect was reported very recently by Hiraoka et al. , 
who produced a monodomain sample of a SmC elastomer 
forming spontaneously from a Srm4 phase upon cooling 
and carried out experiments demonstrating its sponta- 
neous and reversible deformation in a heating an cooling 
process. 

A monodomain SmA elastomer crosslinked in the SmA 
phase is effectively a uniaxial solid with -Dooh symmetry, 
at least for small deformations. For lager deformations, 
however, Srm4 elastomers can show unconventional ef- 
fects, as do Srm4* elastomers in external electric fields. 
These unconventional effects in SmA elastomers are out- 
side the scope of this paper and will be addressed in a 
separate paper [36[ . 

The elastic properties of a SmC elastomer depend on 
whether it was crosslinked in the SmA or SmC phase. If 
it is prepared by crosslinking in a SmC phase, reached 
either by stretching or by applying an external electric 
field, it is a conventional biaxial solid with Cih symmetry. 
If, however, the SmC phase develops spontaneously upon 
cooling from a uniaxial SmA (as in the work of Hiraoka 
et al.) then the underlying phase transition from Dqo/j to 
Cih symmetry (see Fig- HJ) spontaneously breaks the con- 
tinuous rotational symmetry in the smectic planes. As a 
consequence of the Goldstone theorem that requires any 
phase with a spontaneously broken continuous symme- 
try to have modes whose energy vanishes with wavenum- 
bcr, like monodomain nematic elastomers [1, 0, H, @, 0] j 
these SmC elastomers are predicted to exhibit soft elas- 
ticity characterized by the vanishing of a certain elastic 
modulus and the associated absence of restoring forces 
to strains along specific symmetry directions [37| . 

Though biaxial phases are notoriously hard to find in 
nature, it is possible, at least in principle, that a biax- 
ial SmA phase spontaneously forms upon cooling from 
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FIG. 1: Sample distortion and rotations of the Frank director 
n, the uniaxial anisotropy axis e = e 2 , and the layer normal 
N in a transition from a (a) SmA to (b) a biaxial and (c) a 
sheared SmC elastomer. In part (c) we chose the geometry, 
i.e., the coordinate system in target space, so that the smectic 
layers do not rotate. 



a SmA elastomer. In what follows, we will often simply 
refer to biaxial Srm4 elastomers as biaxial smectic elas- 
tomers or biaxial elastomers. In contrast to the afore- 
mentioned phase transition to a soft SmC elastomer, the 
transition to the biaxial SmA phase involves no net tilt 
of the mesogens and takes the system to D^k instead of 
C2/1 symmetry (see Fig. [I]). Nontheless, this transition 
also breaks the rotational invariance in the smectic lay- 
ers spontaneously and thus the emerging biaxial phase is 
soft [37| similar to ideal nematic and SmC elastomers. 

In this paper we study the phase transitions from SmA 
elastomers to biaxial and SmC elastomers and the elas- 
ticity of the emergent phases. As briefly presented in 
Ref. [37]], we set up three different but related models. 
Our first two models involve only elastic degrees of free- 
dom, i.e., they involve exclusively the usual strain tensor. 
Our third model also includes the Frank director speci- 
fying the direction of local molecular order. Each of the 
models is analyzed within mean-field theory, revealing as 
the primary result the soft elasticity of biaxial and SmC 
elastomers [H, |H, [13] alluded to in the forgoing para- 
graphs. 

The outline of our paper is as follows: Section HT1 briefly 
reviews the Lagrange formulation of elasticity theory in 
the context of uniaxial elastomers to establish notation 
and to provide a starting point for our models to fol- 
low. Section [IIII presents our strain-only theory for biax- 
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ial elastomers. We study the transition from the SmA to 
the biaxial state and the elastic properties of the latter. 
We derive the elastic energy density of the biaxial state 
and discuss its softness with respect to certain shears and 
extensional strains. Section llVI contains our strain-only 
theory for SmC elastomers. We investigate the Sim4-to- 
SmC transition and calculate the elastic energy density of 
the SmC phase. Different manifestations of the softness 
of SmC elastomers are pointed out. Section [V] formu- 
lates our theory for SmC elastomers with strain, Frank 
director, and smectic layers. Using the polar decompo- 
sition theorem, we derive transformations between vec- 
tors that transform according to operations on reference- 
space positions x of the undistorted medium and those 
that transform according to operations on target-space 
positions R(x) of the distorted medium, and we formu- 
late the elastic energy for coupled director and strain in 
terms of nonlinear-strain and director fields that trans- 
form under reference-state operations only. In this ap- 
proach, phase transitions can be studied without spec- 
ifying actual orientation in space. We develop a full 
model free energy that includes contributions from the 
crosslinked network, smectic layer compression, and cou- 
pling between the Frank director and the smectic layer 
normal. Then, as above, we study the phase transition 
from the SmA to the SmC phase and the elastic energy 
density of the emergent phase. In addition, we discuss 
the general form of soft deformations and strains in SmC 
elastomers based on rotational invariance in the smectic 
planes and we elaborate on softness under extensional 
strains. Concluding remarks are given in Sec. IVI1 There 
are in total 4 appendixes which contain technical details 
or arguments that lie somewhat aside the line of thought 
of the main text. 



II. LAGRANGIAN DESCRIPTION OF 
UNIAXIAL ELASTOMERS 

As argued above, ideal SmA elastomers are macroscop- 
ically simply uniaxial rubbers, but with nonlinear prop- 
erties that distinguish them from simple uniaxial solids. 
We will employ the usual Langrangian formalism (38l . [39j 
of elasticity theory. Here, we briefly review key elements 
of this formalism in the context of uniaxial elastomers 
to establish notation and to provide some background 
information. 

In the Langrangian formalism mass points in an undis- 
torted medium (or body), which we take as the reference 
space, are labelled by vectors x. Mass points of the dis- 
torted medium are at positions 

R(x)=x + u(x), (2.1) 

that constitute what we call the target space. Both ref- 
erence space points x and target space points R(x) exist 
in the same physical Euclidean space £ where measure- 
ments are made. Thus, R(x) is a mapping from £ to £. 
Both x and R(x) can be decomposed into components 



along the standard orthonormal basis = x,y, z} of 
£: 

x = x t a l , R(x) = J Ri(x)a i . (2.2) 

Here and in what follows, we use the summation con- 
vention on repeated indices unless we indicate otherwise. 
We will also use the convention that indices from the 
middle of the alphabet run over all space coordinates, 
i, j, k = x,y, z. We choose our coordinate system so that 
the z-axis is along the uniaxial direction of the initial 
reference material. Indices from the beginning of the al- 
phabet, a, 6, c, run over x and y only, i.e., over directions 
perpendicular to the anisotropy axis. 

Though reference- and target-space vectors both ex- 
ist in £, they transform under distinct and indepen- 
dent transformation operations. Let Q_ R and O t de- 
note, respectively, matrices describing transformations 
(which we will take mostly to be rotations but which 
could include reflections and inversions as well) in the 
reference and target spaces; then under these transfor- 
mations, x — > x' = O fi x and R(x) — > R'(x) = O r R(x), 
or in terms of components relative to the a-basis 

i^(x) = Or^Rjix), (2.3a) 
x'i = 0_r,v, Xj . (2.3b) 

Unless otherwise specified, we will view O and 0_ T 
as operators that rotate vectors rather than coordinate 
systems. 

Elastic energies are invariant under arbitrary rigid ro- 
tations and translations in the target space and under 
symmetry operations of the references space of the form 
x — > x' = Q^x + b, where b is a constant vector and 
O is a matrix associated with some symmetry element 
of the reference space. Thus, elastic energies are invariant 
under transformations of the form 

R(x)^R'(x')=£ T R(CZ J , 1 x + b) + X, (2.4) 

where is an arbitrary target-space rotation matrix 
and X is a constant displacement vector. In what follows, 
we will generally ignore the displacements X and b. The 
use of CU 1 in Eq. (|2.4j) rather than O r is a matter of 

convention [40]. With the choice O" 1 , when Q_ T = 5 

where 6 is the unit matrix, the mapping R' (x) takes the 
point O r x to the same point in £ as the mapping R(x) 
takes the point x. 

We will usually represent the the reference-space points 
x in terms of their coordinates relative to the basis {a^}. 
We will, however, find it useful on occasion to consider 
orthonormal bases locked to the reference medium and to 
represent reference-space vectors relative to them. The 
initial basis {ej|i = x,y,z} is identical to the basis 
{sii\i = x,y,z}, and x = Xi&i = x^. e z is thus a 
vector along the uniaxial axis of the undistorted body 
[HI]. Under rotations of the body basis, x = x'^ — a^ei, 
where 

ej = R>ij ej (2.5) 
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and x\ = On^jXj. As we shall discuss in more detail 
in Sec. IV Al associated with each reference-space vec- 
tor, there is a target-space vector of the same length. 
Thus, associated with the reference basis {e^}, there is 
a target basis {e^}. In particular there is a target-space 
anisotropy direction e z associated with e z . 

Distortions of the reference medium are described by 
the Cauchy deformation tensor A with components 

Ay = dR t /dxj ee djRi . (2.6) 

It transforms under the operations of Eq. (|2.3|) according 
to 

Afx)-^^')^ 1 (2.7) 

i.e., the right subscript transforms under target-space 
rules and the left under reference-space rules. Usually, 
Lagrangian elastic energies are expressed in terms of 
the Cauchy-Saint-Venant (38l . l42l | nonlinear strain tensor 
u = (g — S)/2, where 

g = A T A (2.8) 
is the metric tensor. The components of u are 

«.jW4( A ^-^) ( 2 - 9a ) 
= | (diUj + djUi + d l u k d j u k ) . (2.9b) 

The strain a is a reference-space tensor: it is invari- 
ant under transformations O t in the target space, but 
it transforms like a tensor under reference-space trans- 
formations: 

i(x)^2 fl ^(x')^ 1 . (2.10) 

This expression applies both to physical transformations 
of reference-space vectors under x\ = 0~j^^Xj or un- 
der changes of basis described by Eq. (|2.5p under which 

To discuss incompressible materials, such as most elas- 
tomers, it can be more appropriate to use variables other 
than the strain tensor to account for deformations that 
are not pure shear. In the case of uniaxial elastomers, 
such variables are the relative change of the system vol- 
ume V, 

V ee SV/V = [det A T A] 1/2 - 1 = [(det(l + 2^)] 1 / 2 - 1 , 

(2.11a) 

and the relative change of separation of mass points 
whose separation vector in the reference state is along 
the z axis: 

r, z ee 8L Z /L Z = lA^A^I 1 / 2 - 1 = (1 + 2u zz ) 1 ' 2 - 1 . 

(2.11b) 

Using these variables and, where appropriate, the ele- 
ments of the strain tensor, the elastic free energy den- 
sity of a uniaxial elastomer to harmonic order can be 



expressed as 

Ami = \ Ci rf z + C 2 VzV +h°3V 2 

+ C 4 u 2 ab + C 5 u 2 az , (2.12) 

where 

u a b = u ab - \5 ab u cc (2-13) 

is the two-dimensional symmetric, traceless strain tensor 
with two-independent components that can be expressed 
as u = ui(e x e x — e y e y ) + U2(e x e y + e y e x ). The elastic 
constant C\ describes dilation or compression along z and 
C3 describes expansion or compression of the bulk vol- 
ume. C*2 couples these two types of deformations. C4 and 
C5 respectively describe shears in the plane perpendicu- 
lar to the anisotropy axis and in the planes containing it. 
With the variables used in Eq. (|2.12[) it is evident that 
the incompressible limit corresponds to C3 — > 00. 

If one approximates r\ and r\ z by the respective leading 
terms in the strains one is left with 

»? = u« + 0(«? i ), (2.14a) 
Vz = u zz + 0(u 2 zz ) . (2.14b) 

and the elastic energy density (|2.12[) reduces to the more 
standard-type expression 

/urn = I C\ u 2 zz + C 2 u zz u u + |C 3 u\ 

+ C 4 u 2 ab + C 5 u 2 az . (2.15) 

Our models to be presented in the following are in spirit 
Landau expansions in powers of Uij or in powers of «y 
and the Frank director, respectively. Hence, for our pur- 
poses, the approximations in Eqs. (|2.14[) will be sufficient, 
and we can use Eq. (|2.15p as a starting point for the 
construction of our models. As we show in Appenix [Al 
using Eq. (|2.15|) , instead of the more general elastic en- 
ergy density Eq. (|2.12p . leaves our results qualitatively 
unchanged, though the more general theory is needed for 
a correct description of the incompressible limit. 



III. BIAXIAL SMECTIC-A ELASTOMERS - 
STRAIN-ONLY THEORY 

In our first theory we consider the case that the shear 
modulus C4 becomes negative as it will in response to 
biaxial ordering of the constituent mesogens of a uniaxial 
SmA elastomer. 



A. Phase transition from uniaxial to biaxial 
elastomers 

If C4 becomes negative, order of the shear strain u ab 
sets in, and higher-order terms featuring u ab have to be 
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added to Eq. (|2 . 1 5[) which leads to the model elastic en- 
ergy density 

ful = Ami + A! u zz u 2 ab + A 2 uuu 2 ab + B (u 2 ab ) 2 , (3.1) 

where we have dropped qualitatively inconsequential 
higher order terms. For the analysis that follows it is 
useful to regroup the terms in by completing the 
squares in \C\u 2 zz + Aiu zz u 2 b , etc., and to reexpress it 
as a sum of two terms, 



= ,(1,1) ,(1,2) 
J uni •> uni * J uni ' 



(3.2) 



where 



& = 3 C rf* + C * v ** v ™ + 3 C * vl + C 5 u 2 az (3.3a) 
f£f =C i u 2 ab + B R (u 2 ab ) 2 



2\2 



2C7 1 (u;+«|)+4B fl (uJ+«5) 



(3.3b) 
(3.3c) 



(1 2) 

The energy /^i ; is clearly identical to the energy of an 
model with a two component vector (ux,U2). Here, 
we have introduced the composite strains 



Vzz 
Vii 



U zz 
Uii - 



(3.4a) 
(3.4b) 

where a and are combinations of the coefficients in fii, , 

' •> uni ' 



b i 



This equation of state is solved by a u ab of the form 



C a Cb 



k&ab 



(3.9) 



where c is any unit vector in the xy-plane and where S 
is a scalar order parameter that takes on the values 



S 







±^/-C 4 /B R 



for C* 4 > , 
for C 4 < . 



(3.10) 



For simplicity, we choose our coordinate system so that 
the x-axis is along c. Exploiting definition (|2. 13|) and 
Eq. (|3.9p and taking c = e x , we find that the equilibrium 
strain tensor m° of the new state for C4 < is diagonal 
with diagonal-elements 



-S+-(p-a)S 2 , (3.11a) 
~S + |09 - a)S 2 , (3.11b) 
\aS 2 . (3.11c) 



,0 



Thus, the new state is biaxial with Dih symmetry. 

The strain u° provides a complete description of the 
macroscopic equilibrium state after the phase transition 
to the biaxial state, but it provides no information about 
a sample's specific orientation in space. The latter infor- 
mation is contained in the Cauchy deformation tensor 







C\ C3 



± ( CzA x - C 2 A 2 

C 2 \ CiA 2 - C 2 A X 



(3.5) 



The subscript R in Eq. (|3.3bj) indicates that the elastic 
constant B is renormalizcd by the ordering of u a b'- 



Br 



B 



1„,2 
2 



a 2 d -al3C 2 - ±f3 2 Q 



(3.6) 



Note that the coefficient (3 vanishes in the limit C3 — > 00. 
a and Br, on the other hand, remain non-zero. We will 
assume that Br remains positive. If it did not, we would 
have to add higher-order terms in u a b, and the transition 
to the biaxial phase would be first order. 

Now we determine the possible equilibrium states u° 

of our model by minimizing As is evident from 

Eq. (|3.3aj) , is minimized for a given equilibrium value 
u° ab of u ab when 





zz 
.0 



as well as 



a ( u °ab) 2 , 

>=P(Uab)\ 

u°=0. 



The equilibrium value of u a b minimizes /, 
tcrmined by the equation of state 



(1,2) 



C i u ab + 2B R u ab {u° cc ,f 



0. 



(3.7a) 
(3.7b) 

(3.7c) 
and is de- 

(3.8) 



A?- = dR° i /dx j , 
which is related to u° via 



4° = \ (A 0T A° 



I) 



(3.12) 



(3.13) 



Note, that the equilibrium deformation tensor A is not 
uniquely determined by vP since rotations in the target 

space change A but do not change u°. Because u° is 
diagonal, it is natural in the present case not to rotate 
the strain after the transition. Then, A is also diagonal 
with diagonal elements given by 



A u 



y/l + S +\{l3-a)S 2 , 



yy 
A 



A° =Ji-S+H(3-a)S 2 , 



Vl + aS 2 . 



(3.14a) 
(3.14b) 
(3.14c) 



It is worth noting that the limit C3 —> 00 in which (3 but 
not a becomes infinite does not yield the incompressibil- 
ity condition det A = 1. This is because in our model, 
C 3 multiplies u\ and not (rj ~ l) 2 [See Eq. (|2~12} ]. 

The emergent anisotropy of the new state in the xy- 
plane can be characterized by the anisotropy ratio 



A 

XX 

A 

yy 



(3.15) 
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Having the equilibrium deformation tensor and the 
anisotropy ratio we can express the scalar order param- 
eter S as 

S = \ [(AL) 2 - (A°,) 2 ] = \ (A° y ) 2 (r ± 1) . (3.16) 

In other words, S is a direct measure for the spontaneous 
anisotropy in the xy-plane. 



B. Elasticity of the biaxial phase 



To determine the elastic properties of the new state, 
we expand in powers of 



8u = u — 



(3.17) 



Since the equilibrium values of v zz , va and u az are zero, 
the expansion of /ii- is trivial, 



<5/uni = \ Cx{5v zz f + C 2 8v„6v« + iC 3 {8vuf 



f (M) _ l 

+ C 5 {8u az } 2 , 
where, up to linear order in duij, 



(3.18) 



5v zz = 8u zz - aS(8u xx - 8u vy ) , 



(3.19a) 

8vu = 8u zz + (1 - [3S)8u xx + (1 + f3S)8u yy , (3.19b) 
8u az = u az . (3.19c) 

As discussed after Eq. (|3.3b|) . the structure of is 
identical to that of an xy model, which has no restoring 
force perpendicular to the direction of spontaneous order, 
which we take to be along the direction. Thus with 
order producing a nonvanishing u\ — (u xx —Uy y )/2, there 
is no restoring force for u 2 — u xy , and 



<*/uni = B R S 2 (Su xx - 8u yy ) 



(3.20) 



Thus Sf^l does not depend on 8u xy to harmonic order, 
and we can conclude already at this stage that the system 
is soft with respect to shears in the xy plane of the orig- 
inal reference material. Merging Eqs. (|3.18|) and (|3.20|) 
and after expressing a, (3 and Br in terms of the original 
elastic constants, we obtain 

Sf {1) = | [Ci + 2C 2 + C 3 ] (8u zz ) 2 + C 5 {8u az f 
+ [d + C 3 + (A 1 + A 2 )S] Su zz Su xx 
+ [C 2 + C 3 - (^i + A 2 )S] Su zz Su yy 



\[C 3 + 2A 2 S + 2BS 2 ] (Su xx ) 
\[C 3 -2A 2 S + 2BS 2 ] (Su yy ) 
[C 3 - 2BS 2 ] 8u xx 8u. 



(3.21) 



after some algebra. 

The strain Su describes distortions relative to the new 
biaxial reference state measured in the coordinates of the 



original uniaxial state. It is customary and more intu- 
itive, however, to express the elastic energy in terms of a 
strain 



«' = (A 0T )- 1 ( 5u(A )- 1 



measured in the coordinates x\ = Xi + i 
new biaxial state. In terms of u' , Sf^ becomes 



(3.22) 
A^Xj of the 



D 2 



c z 



( u zz) "b 2 C X zxz ( u xz ) + 2 C y zyz i u yz ) 



+ Czzxx U ZZ U XX + C zzyy U ZZ U + „ C-. 



2 Cyyyy (u ) 2 + C xxyy U XX U 



yy "zz^yy ~ 2 ^ x %xx 
VV u, xx U yy 



{< x f 

(3.23) 



Our results for the elastic constants Cijki are listed in 
appendix IB H These elastic constants depend on the 
original elastic constants featured in Eq. (|3.1[) and the 
order parameter S and one retrieves the uniaxial elastic 
energy density (|2.15[) for S — > 0. 

Equation (|3.21[) highlights a problem with approximat- 
ing -q = 8V/V with its linearized form uu. In Eq. (|2.12p . 
the limit C 3 — > oo enforces incompressibility, i.e., no 
volume change, even if there is a phase transition. In 
Eq. (|2.15|) . on the other hand, this limit only keeps 
Uu — 0, and uu does not measure a volume change rel- 
ative to a state whose shape has been changed because 
of a phase transition. This is easily seen by noting that 
Sua = A° zz u' zz + A- xx u' xx + A yy u' vv \s not proportional to 
u' u . Thus, the limit C3 ^ oo does not enforce 8V/V = 
in the biaxial phase. If the full nonlinear theory of Eq. 
(|2~T2|) is used, C 3 multiplies (SV/V) 2 in both the uni- 
axial and biaxial phases as we will show in App. [AJ In 
what follows, we will continue to use free energies that 
are harmonic in non-ordering nonlinear strains because 
they give rise to far less algebraic complexity than do the 
more complete theories. The important feature of soft 
elasticity and other physical quantities are not sensitive 
to which theory we use. If detailed treatment of incom- 
pressibility is important, the more complete theory can 
always be used. 

Because there was no 8u xy term in the expansion of 
fiiii there is no term 



C xyxy {u xy ) 



(3.24) 



as there would be in conventional orthorhombic systems, 
because the elastic constant C xyxy is zero. Thus, to lin- 
ear order in the strain, there is no restoring force to xy- 
stresses 



43] 



'xy 



a fsoft 



(3.25) 



i.e., to opposing forces along ±e x applied to opposite 
surfaces with normal along zte^ or opposing forces along 
±e y applied to opposite surfaces with normal along ±e x , 
see Fig. 

Note that there is an interesting parallel between 
the biaxial smectics considered here and biaxial nemat- 
ics formed spontaneously from an isotropic elastomer. 
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dS 



f II -< 



-> f II e. 



|| — e„ 



dS || — 




dS II e r 



to systems described by a Landau expansion of the free 
energy. First, we will consider global rotations in the 
xy-plane of the reference space of our biaxial elastomers. 
On the one hand this will prepare the ground for under- 
standing the softness of extensional strains, and on the 
other hand, it will allow us to understand the vanishing 
of the elastic constant C xyxy from a somewhat different 
perspective. Then we calculate the energy cost of soft ex- 
tensional strains and finally we comment on experimental 
implications. 

We saw in Sec. IIII Al that the direction of the sponta- 
neous anisotropy xy-plane, or in other words, the direc- 
tion of the c-director c, is arbitrary. Thus, equilibrium 
states characterized, respectively, by vP and O z «° 0~ , 
where O „ describes an arbitrary rotation in the refer- 

=R,z J 

ence space about the z axis, must have the same energy. 
With 



FIG. 2: Soft shears in biaxial elastomers. The force f ex- 
erted across an infinitesimal surface element dS is /; = aijdSj. 
There are no restoring forces for (a) external forces along ±e x 
applied to opposing surfaces with normal along ±e s , and (b) 
for external forces along ±e y applied to surfaces with normal 
along ±e x . 



Warner and Kutter [4J] predicted that these biaxial ne- 
matics have many soft modes, and one of these is identi- 
cal to the soft mode of biaxial smectics discussed above. 
Finally, we observe that a biaxial smectic has the same 
point-group symmetry D2h as an orthorhombic crystal. 
However, because it is formed via spontaneous symme- 
try breaking from a state with -Dooh, unlike an equilib- 
rium orthorhombic crystal, it exhibits soft elasticity in 
the icy-plane. An orthorhombic state can be reached via 
a symmetry-breaking transition from a tetragonal state, 
which exhibits square symmetry in the xy plane. Rather 
than exhibiting the soft elasticity discussed above such an 
orthorhombic system exhibits martinsitic elasticity (4f| 
in which domains of different orientation are produced in 
response to stress perpendicular to the stretch direction 
in the xy plane. 



C. Rotational invariance and soft extensional 
strains 

In addition to the softness under u' xy shears, the new 
phase is, provided that the experimental boundary con- 
dition are right, soft with respect to certain extensional 
strains. In the following paragraphs, we will discuss this 
softness in some detail. The mechanism at work here 
is intimately related to a mechanism leading to a sim- 
ilar softness in nematic elastomers, and thus our rea- 
soning will follow closely well known arguments for ne- 
matic elastomers 0, 0|- The form of these soft strains 
depends only on symmetry and the nature of the broken- 
symmetry phase; it is not restricted to small strains or 



a R , z 




(3.26) 



describing a counterclockwise rotation of vectors in the 
reference space through ■& about the z-axis, the strain 

»'(*) - (A° T ) 1 \Q. R z M° - M°] (A°) 1 (3.27a) 
, , -rj/(l-cos2z9) r~ 1/2 sin2i? 



,-1/2 



sin 2$ 







1 - cos 2tf 




(3.27b) 



must not cost any elastic energy. Our elastic energy den- 
sity (|3.23[) contains only second-order terms in u--, and 
hence it can at best be invariant with respect to infinites- 
imal rotations (46j. However, even for infinitesimal the 
strain w'(i9) has nonzero components, namely 



1 



(3.28) 



Thus, as it does in Eq. (|3.23[) . the elastic constant C xyxy 
of the term (|3.24p must indeed vanish, and this vanishing 
can be understood as a result of the spontaneous sym- 
metry breaking in the xy-planc. 

The existence of zero-energy strains that reproduce ro- 
tations in the reference space has consequences reaching 
further then just the softness with respect to shear strains 
u' xy , viz. depending on the experimental boundary con- 
ditions, extensional strains u' xx and u' yy can also be soft. 
If the boundary conditions are such that no relaxation 
of strains is allowed, then strains u' xx and u' yy will cost 
an elastic energy proportional to {u' xx ) 2 and (u yy ) 2 , re- 
spectively. If, however, strain relaxation is allowed and 
one imposes for example u' yy with the right sign, then 
u' xx and u' can relax under the right circumstances to 
produce the zero-energy strain of Eq. (13.27p . i.e., to make 
u' yy a soft deformation. 
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To discuss this in more detail, let us assume for con- 
creteness that the anisotropy in the xy-plane is positive, 
rj_ > 1. Then it follows from Eq. (|3.27p that u' yy is pos- 
itive and that u' xx is negative for soft strains and that 
we consequently can only have soft elasticity for u' yy > 
and u' xx < 0. Let us consider here as an example a strain 
with u' yy > 0, i.e., a stretch of the sample along y. Com- 
parison with Eq. (|3.27|) shows that, if strain relaxation is 



allowed and u' . 



and u' relax to 



-r_,_ Uy y - 



~xy 



±1 



'< y (rs 



1 - 2u>. 



2r ± 



(3.29a) 
(3.29b) 



then u'yy is converted into a zero-energy rotation through 



an angle 



1 

±— sin 
2 



r± 



-l\j 2u 'yy( r ^- l - 2u 'yy) 



(3.30) 



Thus, u' costs no elastic energy as long as < u' yy < 
(r± - l)/2. 

When u' yy is increased from zero, $ increases from zero 
until it reaches ir/2 at u' yy = (r± — l)/2 [and u' xx — 
(rj 1 — l)/2, u' xy = 0] at which point, the deformation 
tensor defined by A^ - = dRi/dx', is 



-1/2 
_L 











1/2 











1 



(3.31) 



which leads via Ay- = A^A'^ to an overall deformation 







A = 



A 

yy 

Ay y i/rZ 
A° 



(3.32) 



relative the original uniaxial state. Thus, at $ = tt/2, Ay- 
is identical to A°- except with A^ = ^/r±A yy replaced 
by A yy and A yy replaced by A^, i.e., the x and y axes 
have been interchanged in going from •& = to i? = tt/2. 
In the process, c rotates from being parallel the x-axis to 
being parallel to the y-axis • 

For u' yy > (r± — l)/2, there is no real solution for 
i9, and a further increase in u', measured by Su' yy = 
u'yy — (r± — l)/2, which stretches the system along the 
space-fixed y-axis, cost energy. Since the anisotropy axis 
is now along the y- rather than the x-axis, this stretching- 
costs the same energy as it would have cost to stretch 
the original system with anisotropy axis along the space- 
fixed x-axis by the same amount. The yy-component of 
the strain relative to the state with $ = tt/2 is &u' yy = 
(A' 0yy )~ 2 6u' yy = rj_ Su'y. Thus, because the x- and y- 
axes have been interchanged, the free energy as a function 
of u' yy is 



/•soft 

Jn 2h 



for Su' yy < , 



\Y x (5u') 2 for Su' >0. 



(3.33) 




FIG. 3: Schematic plot (arbitrary units) of the engineering 
stress <r™ g versus the deformation A' yy for a soft biaxial elas- 
tomer with an equilibrium order in the xy-plane along x and 
positive anisotropy, r± > 0. Up to a critical deformation 
A' 0yy — y[r~Z the sample responds to the deformation merely 
by rotating c and consequently the stress is zero. Above 
Ao TO the sample stretches along the new direction y of the 
c-director and <r™ s grows linearly for small A' yy — A' 



Oyy 



where 



Yr 



Cxxxx 



^yyyy^zzxx 



C 'zzzzC ' xxyy 2C X xyyC X xzzCy 



C C — 

^yyyy^zzzz ^yyzz 



(3.34) 



is the is Young's modulus for stretching along the 
anisotropy axis in the xy-plane (originally along x). 

Equation (|3.33[) has tangible implications for the 
stress-strain behavior of soft biaxial elastomers. The 
stress that is usually measured in experiments is the engi- 
neering stress, i. e., the force per unit area of the reference 
state. For the extensional strain under discussion here, 
the engineering stress [43[ is to leading order 



fi £ soft 

dA 



mi 





YxK 



for A'yy < ^JTZ 



Qyy ^'yy fOr Ayy 



(3.35) 



For 



cng 

yy 



a 



i' 

yy 

*Y X (A> 



mi 



.. .. KyyiKv ~ A 'oyy)> and 
Aq ). Figure [3] depicts the dependence 



A' 



6u 'yy 



of (Tyy S on A' yy . o m A^ y — up to a, critical deforma- 
tion A' 0yy = ^/r± the stress is zero. Above the critical 
deformation, a™ s grows linearly from zero. 



IV. SMECTIC-C ELASTOMERS 
STRAIN-ONLY THEORY 



In this section we use the strain-only theory to study 
the phase transition from a uniaxial SmA elastomer to 
a SmC elastomer when C5 becomes negative in response 
to a SmC-ordering of the mesogenic component. 
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A. Phase transition from uniaxial to smectic-C 
elastomers 

When C5 is driven negative, the uniaxial state becomes 
unstable to shear in the planes containing the anisotropy 
axis, and the uniaxial energy (|2. 15[) must be augmented 
with higher-order terms involving u az to stabilize the sys- 
tem: 

/uni = /uni + Di U zz U 2 az + D 2 Unu\ z + D 3 U ab U az U bz 

(4.1) 



+ E(u 2 az ) 2 , 



where we omit all unimportant symmetry-compatible 
higher-order terms and we use u zz and Ua rather than 
r\ z and 77. To study the ordered phase of this free energy 
when C5 < 0, we proceed in much the same way as we 

did for the biaxial state of f^. Using Eq. (12.15|) for / un i, 

(2) 

we complete squares to write / uni as the sum 

,(2) = ,.(2,1) f (2,2) „v 

J uni •> uni ' J uni ' V / 

of the two terms 

finP = w lz + C 2 w u w zz + \ C 3 w% + C 4 w 2 ab , 

(4.3a) 

(4.3b) 



/<!f> = C 5 u 2 az + E R (u 2 az ) 2 



where E Rl which we assume to be positive, is a renor- 
malized version of E, 

E R = E- \a 2 & -arC 2 - \r 2 C 3 - \lo 2 C 4 , (4.4) 

and where 



w zz = u zz - ou 2 az , (4.5a) 
W U = Uii - Tu 2 az , (4.5b) 
Wab =u ab -uj (u az u bz - \8 ab u 2 cz ) . (4.5c) 

The coefficients in Eqs. (14.4)) and (|4. 5[) are given by 



G 



— 1 



C3D1 — C2-D2 

d C 3 - C 2 \CiD 2 - C 2 D 1 



and 



D3 
2C 4 



(4.6) 



(4.7) 



In the limit C3 — » 00 the coefficient t vanishes whereas 
tr, to, and Er remain finite. 

The equilibrium value u® z of u az is determined by min- 

..(22). 
imizing / un - , which has xy symmetry. Provided that we 

choose our coordinate system so that its ir-axis is along 

the direction of ordering, the corresponding equation of 

state, 



C 5 u° az + 2E R u az (u° bz ) 2 = Q, 



(4.8) 



leads to 



and 



u° z = for C 5 > and C 5 < (4.9) 



S — u xz 



for C 5 > . 



±y/-C 5 /(2E R ) for C 5 <0. 



(4.10) 



From Eqs. (|4.3a[) and (|4. 5[) the other components of u° 
follow as 



U xx = 2~ {t + uj — a) S 2 , 
u vv = \ ( T ~ w ~ a ) 52 ' 



u 



vv 




aS 2 



(4.11a) 
(4.11b) 
(4.11c) 



We learn that, unlike the biaxial case, vP is not diagonal; 
it has nonvanishing xz and zx components, which implies 
C 2 h rather than D 2 h symmetry. 

As already pointed out in Sec. lIIIl an equilibrium strain 
tensor determines the corresponding equilibrium defor- 
mation tensor only up to global rotations in the target 
space. Here we choose our coordinate system in the tar- 
get space so that the transition from the D 2oo to the 
C 2 h state corresponds to a simple shear as shown in 
Fig. [U i.e., we choose the target space coordinates so 
that tan^ = A° z /A° z is nonzero but A^ = 0. Then the 
only nonzero components of A° are 



A° xx = y/l + 2u° xx = ^/1 + (t + uj-<j)S 2 , (4.12a) 
A° y = yjl + 2u° yy = ^l + (r-io-a)S 2 , (4.12b) 

A°, = 2 ^- ■ (4.12c) 



9 a xz 

A° 



V^l + (t + w - <t) S 2 



y/l + 2u\ z - (A° xz y 



=\A +2 (a "ifr + i r+ 7; )aS4 ' ^ 

y 1 + (r + lu - a) S z 

B. Elasticity of the smectic-C phase 

Next we expand the elastic free energy density about 
the equilibrium state. The expansion of /^f [Eq. (|4.3ap ] 
is particularly simple and leads to 

SfunP = I Ci (Sw zz ) 2 + C 2 Sw u Sw zz + |C 3 (Swu) 2 

+ d (Sw ab ) 2 (4.13) 

with 



5w zz = 6u zz - 2aS 5u xz , 
Sw u = Sua - 2tS 5u xz , 



(4.14a) 
(4.14b) 



Sw xx = -Swyy = \ (Su xx - Suyy) - lu S Su xz , (4.14c) 

8w xy = 5Wy X = Su xy — COS SU yz . (4.14d) 
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Expanding /£f ) [Eq. (|Ob|) ], we find that 

<5/f„'i 2) = ±E R S 2 {8u xz ) 2 (4.15) 

is independent of 8u yz , and we might naively expect the 
system to exhibit softness with respect to u yz . This, how- 
ever, is not the case because f^P depends on 5u yz via 



the relative strain l|4.14d[) . Thus, the softness of the or- 
dered phase with C 2 h symmetry is more subtle than that 
of the biaxial phase with D 2 h symmetry, as we will dis- 
cuss in more detail further below. Assembling the contri- 
butions (|4.13p and (|4. 1 5[) we obtain for the entire elastic 
free energy density to harmonic order 



SfW = 2C A [5u xy - uS6u yz } 2 + i [d + 2C 2 + C 3 ] {5u zz f + AES 2 (5u xz ) 2 + \ [C 3 + C 4 ] {{5u xx f + {Su yy ) 2 } 
+ [C 2 + C' 3 ]Su zz {Su xx + Su yy } + [C 3 - Ci\Su xx 5u yy + (2D 2 + D 3 )S Su xx Su xz + (2D 2 - D 3 )S 5u yy Su xz 
+ 2(D 1 +D 2 )SSu zz Su xz . (4.16) 



Our remaining step in deriving the elastic free energy 
density of SmC elastomers is to change from the variables 
of the old uniaxial state to those of the new state by 
switching from Su to the new strain tensor u' as defined 
in Eq. (|3.22[) . With the equilibrium deformation tensor 
as stated in Eqs. (|4.12p we obtain after some algebra 

fbf h = \ C [cos 9 u' xy + sin 9 u' yz ) 2 + \ C zzzz (u' zz ) 2 

+ 2 C'xzxz i u xz ) + Czzxx U zz U xx ~l~ ^zzyy u zz u yy 
+ 2 C'xxxx ( u xx ) + 2 C yyyy ( u yy ) + Cxxyy u xx u yy 

+ C xxxz u' xx u' xz + C yyxz u' yy u' xz + C zzxz u' zz u' xz , (4.17) 

where the angle 9 and the elastic constants depend on 
the original elastic constants in Eq. (|4.ip and S so that 
one retrieves the uniaxial energy density (|2.15jl for S — > 
0. The explicit results for these quantities, which are 
somewhat lengthy, can be found in Appendix [B] In the 
incompressible limit, the specifics of 9 and the elastic 
constants get modified, however, without changing the 
form of Eq. gT7D. 

Having established the result (|4.17[) . we are now in the 
position to discuss the anticipated softness of the new 
state. Our first observation is that the elastic energy 
density of Eq. (|4.17[) has only 12 (including 9) rather 
than the 13 independent elastic constants of conventional 
monolinic solids [48| . That is because here there are only 
two rather than three independent elastic constants in 
the subspace spanned by u' xy and u' yz . Below, we present 
two derivations of this result. 

In the first derivation, we exploit the fact that 
cos 9u' xy + sin 6u' can be viewed as the dot product of 
the "vectors" v — (u' x y ,u' yz ) and e\ = (cos 9, sin 9). Thus, 
the first term in Eq. (|4.17|) . 

\C [cos9u' xy + sin9u l yz ] 2 = \C{e 1 -v) 2 1 (4.18) 

is independent of e 2 ■ v, where 

e 2 = (-sin0,cos0) (4.19) 

is the vector perpendicular to e% . Thus, distortions of the 
form — sin 9u' + cos 9u yz = e 2 ■ v, i.e, distortions with v 



along e*2, cost no elastic energy. A manifestation of this 
softness is that certain stresses cause no restring force and 
thus lead to large deformations. To find these stresses 
we take the derivative of the elastic energy density (|4. 17|) 
with respect to u' xy and u' yz which tells us that 

— sin 9 a xy + cos 9 a zy — e 2 ■ w — , (4.20) 

where again, cry is the second Piola-Kirchhoff stress ten- 
sor and where we have introduced the "vector" 

w = {(T xy ,a zy ) (4-21) 

in the xz-plane. Equation (|4.20p means that there are no 
restoring forces for external forces in the rrz-plane along 
±e 2 = ±(— sin#, 0, cos#) applied to opposing surfaces 
with normal along ±e y , and there no restoring forces for 
external forces along ±e y applied to surfaces with normal 
along ±e 2 . These kinds of stresses are depicted in Fig. S) 

A second derivation of the softness of the SmC phase 
involves a transformation to a rotated coordinate sys- 
tem, as described in Eq. (|2.5p . in which the first term in 
Eq. (|4.17p is diagonal. The components x", of a reference- 
space vector x' expressed with respect to a rotated basis 
{e'^li' = x',y',z'}, where e^, = R ^ V jhj with 

(cos 9 sin^ \ 
10 (4.22) 
— sin 9 cos 9 J 

describing a counterclockwise rotation of the reference- 
space basis about the y-axis, are simply x" — On y: ij x 1 -. 
The components of the strain matrix expressed in the ro- 
tated basis are u-'y = e'i> • ^ • &'j> = O^inUijO^.^,, 
from which we obtain via Eq. (|2.10p that u x , , = 
cos9u' xy + sm.0u' yz and u y , z , = —sin9u' xy + cos9u' yz . 
Thus, taking 9 to be the angle appearing in Eq. (|4. 17|) 
and dropping the double-prime from the strains, we ob- 
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FIG. 4: Soft distortions of SmC elastomers. There are no restoring forces (a) for external forces in the zz-plane along ±62 
applied to opposing surfaces with normal along ±e a , and (b) for external forces along ±& v applied to surfaces with normal 
along ±62. 



tain 

fc-ih = 2^ x 'v' x 'y' i u x'y') 2^ z ' z ' z ' z ' ( Uz ' z ') 

2 ^x' z'x'z' ij^x' 2') "t" Cz'z'x'x' 'U'z'z' l^x'x' 

H~ C Z l Z l yl y' U z l z'Uylyl ~\~ 1^ C x' X' X 1 X 1 {^x'x'^) 

2 ^v'y'v'y' (^y'y') ^x'x'y'y' ^x' x' '^y' y' 
"T" Cx'x'x'z 1 ^x' x'^x' z' "T" Cy'y'x'z' Uy'y'^x' z' 

+ C z < Z 'x'z' u z < 7j Ux> Z ' ■ (4.23) 

for the elastic energy density of SmC elastomers in the 
rotated coordinates. Note that, in this coordinate sys- 
tem, the elastic constants C X 'y> y 'z' and C y i z i y ' Z ' are zero, 
and the elastic energy does not depend at all on u v i z >. 

Cx'y'x'y' is identical to C and ^y'y'y'y' ~~ ^yyyy ■ The 

remaining new elastic constants are non-vanishing con- 
glomerates of the elastic constants defined via Eq. (|4.17|) 
and sines and cosines of 6. We refrain here from stating 
further specifics to save some space and because calcu- 
lating these specifics is a straightforward exercise. 

The vanishing of C x >y' y 'z' and C y i z > y i z i means that 
SmC elastomers are soft with respect to shears in the 
y'z'-plane. If one can cut a rectangular sample with 
faces perpendicular to the , then there are no restoring 
forces for external forces along ±e z / applied to opposing 
surfaces with normal along ±ey , and for external forces 
along ±ey applied to surfaces with normal along ±e' z ,. 
This softness can be visualized as in Fig. [2] with e x re- 
placed by e z , and e y replaced by e y , = e y . 

Before we go on, let us finally comment on the im- 
pact of the softness on the phonon spectrum of SmC 
elastomers. A comprehensive discussion of the phonon 
spectrum of course requires a dynamical theory of smec- 
tic elastomers. This is beyond the scope of this paper 
and makes the topic of a separate publication [^H . Our 
current theory allows us to address static phonons. To 
this end we switch to Fourier space [49[ and re-express 
the elastic energy density (|4.17[) in terms of the Fourier 
transform u(q) with respect to x' of the displacement 
field u'(x') = R'(x') - x'. Neglecting the non-linear part 
of the strains, cf. Eq. (|2.9[) . this yields a Fourier trans- 
formed elastic energy density f(? h harmonic in u and the 



wavevector. With u x — u z = q y = 0, * reduces to 

fbll = 5 & (cos q x + sin 9 q z f u\ . (4.24) 

Hence, the energy cost is zero for phonon displacements 
u = (0,u y ,0) parallel to e y with wavevector q || e z > . 
By similar means one also finds that there is no energy 
cost for u || e z i with q || e. y . This softness of these 
phonons has tangible implications on the dynamics in 
that it leads, e. g., to a vanishing of the corresponding 
sound velocities [341 ] . 

An interesting question that we have left aside so far is, 
of course, whether SmC elastomers are, like nematics or 
biaxial smectics, soft with respect to certain extensional 
strains. We will postpone this question to Sec. [V] until 
after we have developed our strain-and-director model for 
SmC elastomers. This will then allow us to discuss soft 
deformations and strains more comprehensively including 
their impact on the director. 

V. SMECTIC-C ELASTOMERS - THEORY 
WITH STRAIN AND DIRECTOR 

Our second theory for SmC elastomers, to be presented 
in this section, explicitly accounts for the smectic layers 
and for the director n. It generalizes the achiral limit of 
the continuum theory by Terentjev and Warner [26| in a 
formalism that ensures invariance with respect to arbi- 
trary rather than infinitesimal rotations of both the direc- 
tor and mass points. We will see as we move along that 
the properties of the transition to the SmC phase pre- 
dicted by this theory are identical to those of the strain- 
only theory, discussed in the preceding section, in which 
C5 goes to zero. 

A. Reference- and target-space variables and the 
polar decomposition theorem 

In traditional uncrosslinked liquid crystals, there is no 
reference space, and all physical fields like the smectic 
layer-displacement field U , the layer normal N, and the 
Frank director n are defined at real, i.e., target-space 
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points R, and they transform as scalars, vectors, and 
tensors under rotations in the target space. In the La- 
grangian theory of elasticity, fields are defined at refer- 
ence space points x, and they transform into themselves 
under the symmetry operations of that space. To de- 
velop a comprehensive theory of liquid-crystalline elas- 
tomers, it is necessary to combine target-space liquid 
crystalline fields and reference-space elastic variables to 
produce scalars that are invariant under arbitrary rota- 
tions in the target space and under symmetry operations 
of the reference space. This requires that we be able to 
represent vectors and tensors in either space Q. 

To be more specific, let b be a target-space vec- 
tor, which by definition transforms under rotations to 
b' = 0_ T h, and let b be a reference-space vector, which 

transforms to b = O b. Recall that both reference- and 
target-space vectors exist in the same physical Euclidean 
space E. Therefore, there must be a transformation 
that converts a given reference-space vector to a target- 
space vector and vice versa while preserving length. This 
transformation is provided by the deformation matrix A 
and the matrix polar decomposition theorem [50| , which 
states that any non-singular square matrix can be ex- 
pressed as the product of a rotation matrix and a sym- 
metric matrix. If b is a reference-space vector, then Ab is 
a target-space vector that transforms under O but does 
not change under O because under x — * x' = 0„x and 

IX it 

R -» R' = £ t r, b -> b' = Qjo, a -► A' = QtAPb 1 

and Ab -> A'b' = O AO^O b = £ T Ab. The transfor- 
mation b — > Ab, however, does not preserve length. To 
construct a transformation that does, we simply multiply 
A by the square root of the metric tensor to produce 

£ = As~ 1/2 - (5.1) 

This operator clearly satisfies T = OO t = 6 and 
detO = 1, and it is thus a length- preserving rotation 
matrix. Equation (|5.ip . which can be recast in the form 
A = Qg 1 / 2 is simply a restatement of the polar decom- 
position theorem because g is a symmetric matrix. To 

first order in dui/dxj, Oij reduces to the standard ex- 
pression for an infinitesimal local rotation of an elastic 
body through an angle fl = x u, 

= Sij ~ eijk fijH , (5.2) 

where e^k is the Levi-Civita tensor. Equipped with O 

we can convert (or rotate) any reference-space vector b 
to a target space vector b via 

b = Qb (5.3) 

and a target-space vector to a reference space vector via 

b = Q T ■ b . (5.4) 

An alternative interpretation of the relation between b 
and b follows from 

b = b t e { = O.J'.v, = bjtj, (5.5) 



with 

4- r\T -1/2 <5R , , 

t =O^ 3 =g lk ' — , (5.6) 

where we used A^e; = (dRi/ dxk)&i = dR/dx k . The set 
of vectors {tj|i = x,y,z} forms an orthonormal target- 
space basis in the tangent space of the deformed medium 
(recall that dH/dxj is a tangent-space vector). Thus, hi 
represents the components of the target-space vector b 
relative to the orthonormal tangent-space basis defined 
by U. 

We can now apply this formalism to the Frank director 
in smectic elastomers. The familiar director n is a target- 
space vector, which we can represent as 

n={c,n z ), n z = y/l - c 2 a , (5.7) 

where c is the so-called c-director. Contractions of the 
components of n with the stress strain tensor Uij do not 
produce a scalar because rij and transform under dif- 
ferent operators. To create scalar contractions, we can 
convert n to a reference-space vector via n = _1 n, 
where O is defined by Eq. (|5.ip , with 



h=(c,n z ), n z = y/l - cl . (5.8) 

Combinations like n a u a bnb and n a u az n z are now scalars. 
When linearized, these combinations reproduce those in 
the original de Gennes theory [5lJ . Linearized deviations 
of the target-space director from its equilibrium no can 
be expressed as <5n = u x no, where w is a rotation 
angle. Then 5n = c = (u> — £1) x no and, for example, 
c a Ua Z ~^ u az [(tjj — ft) x no] . Since we are interested in 
the transition from the SmA to the SmC phase in which 
n undergoes a rotation through a finite rather than an 
infinitesimal angle relative to its equilibrium in the SmA 
phase, we need to use the full nonlinear representation of 
rotation matrices. 

As we discussed in Sec. [TT1 the reference space can be 
endowed with an orthonormal basis {&;}. This space 
is anisotropic, and we take e z to be along the uniaxial 
anisotropy direction of the SmA material [4l| . Crosslink- 
ing in the SmA phase freezes in an anisotropy direction 
in the elastic network and, therefore, a general prefer- 
ence for the reference-space director to align along e 2 . 
This preference, present in elastomers crosslinked in the 
nematic as well as the Sm^4 phase, is distinct from the 
preference, which we will discuss shortly, for the director 
to adopt a preferred angle relative to the layer normal. 

An important property of O is that it reduces to the 
unit matrix when A is symmetric, i.e., under pure shear 
transformations, target- and reference-space vectors are 
identical. Thus if a reference-space vector is known (cal- 
culated, for example by minimizing a free energy that de- 
pends only on reference-space vectors and tensors), the 
associated target-space vector is obtained by rotating the 
reference-state vector by the operator O, which is the 
same operator that rotates the pure shear configuration 
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(c) (d) 

FIG. 5: Schematic representation of distortions in the xz- 
plane induced by the SmA-to-SmC transition: (a) undistorted 
Srm4 phase, (b) SmC phase with a symmetric deformation 
tensor A, (c) SmC phase with A xz > and A zx = 0, and (d) 
SmC phase with A xz — and A zx > 0. 

to the target-space configuration described by A. Fig- 
ure [5] depicts the effect on an initial reference state unit 
vector of a symmetric shear and then a subsequent rota- 
tion to a final target state. 

B. Development of a model for smectic elastomers 



can conveniently be decomposed as follows: 

/ /uni F /nonlin "F fc "F /coupl- (^•^) 

We discuss each of these terms individually. / un i is the 
uniaxial elastic energy to quadratic order in strain of 
Eq. (|2. 15[) with 5 elastic constants C\,...,C§. Nonlinear 
strain energies are contained in / n0 niin. This energy has 
many terms in general, but we will keep only those terms 
that are relevant to the development of shear strain and 
smectic-C order: 

/nonlin = ~B X U„t& F B 2 (t&) 2 . (5.10) 

To impose the full nonlinear incompressibility constraint 
to order (h^ z ), we should include couplings of uu to u 2 z . 
If we did this, there would be a contribution to the energy 
of the form ^B[uu — 2u 2 a ] 2 where B is the bulk compres- 
sion modulus. This would of course yield a contribution 
to B2 of order B. To treat such a term, we could replace 
uu by 2u 2 a + Sua and express the energy in terms of Sua 
rather than uu. We will instead continue to consider the 
theory defined above without coupling between uu and 
u 2 z . The two theories effectively differ only in the value 
of i?2- The energy associated with the development of 
c-order described with nonzero c a is 

fc = \cl+ V -{ci)\ (5.11) 

and the energy of director-strain coupling is 

/coupi = Ai c 2 a u zz + A 2 c 2 a uu + A 3 c a u ab c b 

+ A 4 c a u az + X 5 u zz u az c a . (5-12) 



Having established the relation between reference- and 
target space vectors, we can now develop a complete phe- 
nomenological energy density / for smectic elastomers in 
terms of reference-space variables only. For briefness, we 
will in the following often refer to energy densities sim- 
ply as energies. We will be content with an expansion 
of / in powers of the strain Uij and the c-director up to 
fourth order. There are three distinct contributions to /: 
(1) the elastic energy f net of the anisotropic network - 
this is the harmonic energy of nematic elastomers 0, [52j 
augmented by some nonlinear terms, (2) the compression 
energy /i aye r of smectic layers, and (3) the energy /tut as- 
sociated with tilt of the nematic director relative to the 
layer normal. The latter two contributions are essentially 
the Chen-Lubensky (CL) model [H| model generalized to 
elastomers. We will assume that smectic layers are locked 
to the elastic network as is the case when the network is 
crosslinked in the smectic phase. The moduli and cou- 
pling constants in / not , /i ay cr, and /tilt contribute to the 
moduli and coupling constants in /. We will identify and 
calculate these contributions in what follows, denoting 
each contribution by the appropriate superscript. Thus, 
for example Cf et is the contribution of C\ from / net . 

The sum of all of the above contributions to the energy 



Here we have retained terms linear in U{j and up to 
quadratic order c a . We also keep the A5 term, which 
is linear in c a and quadratic in u; this contributes to the 
behavior of SmA elastomers under uniaxial extensional 
strain, which we will discuss in a separate paper [3lj | . The 
contributions of / ne t, /layer, and /tut to the various terms 
in /uni, /nonlin, fc, and / coupl are summarized in Table |U 

We now discuss the individual contributions / ne t, 
/ layor , and /tut, beginning with / net . This is the en- 
ergy of a nematic elastomer, which we express in terms 
of reference-state variables only and which we expand 
in powers of Uij and c a . It can be calculated from 
the neo-classical energy developed by Warner and Ter- 
entjev [l| with an explicit volume compression energy 
i_B(det<7 — l) 2 , where B is the compression modulus 

of order 3 x 10 9 Pa, in addition to the entropic energy 
■i/iTrAi o A T £~ 1 , where /x ~ 10 6 Pa is the rubber shear 
modulus and £ and I are the polymer step-length ten- 
sors, respectively, at sample preparation and in a general 
distorted state of the system. The first step-length tensor 
is a reference-space tensor with components 

£oij = + (j? - l^e,-] , (5.13) 
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TABLE I: Contributions to coefficients in / from /net, /layer and ftnt with a = /j,(p — l) 2 /(2p). 





Ci 


c 2 


c 3 


c 4 


C 5 


Si 


#2 


Ai 


A 2 


A 3 


A 4 


A 5 


r 


V 


net 




-A* 


4B -n 




i„b>+i) a 

2^ p 


a 


—a 




2p 


r V 


r V 


-!« 







layer 


4B sm 














6-Bsm 


9 D 


2-Ssm 











4B sm 







tilt 
















4* 


\rt + \v t 











r t 


-Art 


rt 


Wt 



where e = e z specifies the direction of uniaxial anisotropy 
in the references state (denoted by n in WT), £± is the 
step length perpendicular to e and p = £u/£x> with £|| 
the step length parallel to e. £ _1 is a target-space tensor 
with components = ^J 1 ^ + — l)nj7ij]. Fol- 
lowing standard procedures l[, / nGt can be cast as the 
sum of a uniaxial energy the form of / U m + ./nonim + /cou P i- 
Its contributions to the various moduli and coupling con- 
stants are listed in tablelH The coefficients C§ et , AJ et and 
r nct satisfy Cg ct = (A? et ) 2 /(2r not ) as is required for soft 



nematic elastomers [lj. We should logically add a semi- 
soft energy 0, [f| [fjj because we assume that the system 
was crosslinked in the smectic phase. This will turn out 
to be unnecessary because /i aye r and ftnt contribute the 
same kind of semi-soft terms but with greater magnitude. 

To derive both /i aye r and ftnt, we need to discuss in 
more detail the smectic displacement field U and the layer 
normal N. The smectic mass-density- wave amplitude for 
a system with layer spacing d has a phase 



(R) = q Q [R z - U(R)} , 



(5.14) 



where qo = 2n/d. Since there is a one-to-one mapping 
from the reference-space points x to the target-space 
points R(x), we can express <f> as a function of x as 



(x) = gt,[z + «,(x)-tf(R(x))]. 



(5.15) 



We are only considering systems crosslinked in the smec- 
tic phase in which the smectic mass-density wave cannot 
translate freely relative to the reference material, and 
there is a term 



Using the polar decomposition theorem, we can calculate 
the reference-space layer normal 



Ni = 



[9 



-1/21 



-llV2 

\zz 



(5.19) 



With these definitions, we have 



n • N = n • N w 1 



l{ul z f. 



2 ^zzCaUaz 



(5.20) 



In this expression, we have retained the dominant terms 
necessary to describe the SmA-SmC transition and the 
Helfrich-Hurault instabilities 0, [EBJ produced by an 
extensional strain u zz along z. We have not included 
higher-order terms in c a and Uij , which could change the 
numerical values of our results but not their form. 

The preferred spacing between smectic layers depends 
on the orientation of the director relative to the layer nor- 
mals. If n is parallel to N, the preferred spacing is d. If 
n is not parallel to the smectic layer spacing should scale 
approximately as dcos0, where is the angle between 
n and N. A phenomcnological energy that reflects this 
preferences is 



/layer = ~B sm q Q 4 [(n ■ V0) 2 



- 2 B sm 



2£L 



[<r%(n-N) 2 -l 

Uzz + \{c a + u az ) 2 - 2u 



1 12 
az\ ■ 



(5.21) 



/lock-in = I A (U Z - U) 2 



(5.16) 



in the total free-energy density that locks the smectic 
field U to the displacement field u z [27| . In what follows, 
we will take this lock-in as given and set U = u z . This 
has some interesting consequences. The smectic phase is 
now <p = qoz, which implies 



Vi0 = 00/aRi = g o [A- 



(5.17) 



where we introduced the notation that is the ij- 

component of the matrix M a for any matrix M_ and ex- 
ponent a (we retain the notation My = [M]y ). Thus, in 
the target space, the unit layer normal reads 



Ni 



[A- 1 ]. 



(5.18) 



The smectic compression modulus B sm is of order 10 7 Pa 
deep in the smectic phase though it vanishes as the tran- 
sition to the nematic phase is approached. /i ay0 r is the 
generalization to elastomers of the compression energy in 
the Chen-Lubensky [53| model for SmA and SmC phases. 
We could have used the more isotropic compression en- 
ergy proportional to [(V(f>) 2 — q 2 ,] 2 instead. It is the one 
studied by AW [3l|. The advantage of the CL energy 
over the latter energy is that it encodes the tendency for 
layer spacing to decrease when becomes nonzero. If, 
for example, u az = and c a ^ 0, then to minimize /layer, 
u zz — —\(? a < 0. Thus, as expected tilt decreases u zz 
and layer spacing. The CL energy, like the more isotropic 
one, also has built in the physics of the Helfrich-Hurault 
instability (54|, which we will discuss in another paper 
[36| . If c a = and u zz ^ 0, then /i ayor is minimized 
when there is a shear strain, u za — ±yJu zz /2. 
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Finally, we turn to the tilt energy. This is most easily 
expressed in terms of sin 2 = 1 — (n • N) 2 : 

1 „ 1 , 
/tilt = -n sin 6 + -v t sin 9 



r t [(c a + u az ) 2 - u zz u az c a - u zz u 2 az + (u 2 az ) 2 } 

(5.22) 



+ ^v t [(c a + u az ) 2 ] 2 . 

The modulus rt is generally of order but less than B sm . 
However it vanishes as the transition from the SmA to 
the SmC phase in uncrosslinked smectics is approached. 



C. Phase transition from smectic- A to smectic-C 
elastomers 

Having developed a full model for smectic elastomers 
that provides a description of both the SmA and SmC 
phases, we can study the transition from the SmA to 
the SmC phase. In this transition c Q becomes nonzero, 
and because of the coupling between c a and u az , u az 
also becomes nonzero. Alternatively, we could say that 
the angle ~ c a + u az becomes nonzero and drives the 
development of a nonzero u az because of a — u az cou- 
pling. We will use the variables c a and u az to describe 
the SmA-SmC transition. To keep our discussion sim- 
ple, we will focus on the development of c-order and in- 
clude only those terms in the free energy that play an 
important role in this transition. Accordingly, we will 
ignore / non iin, i-e., we set B\ = B% = 0), and we will set 
A5 = 0. Setting these coefficients, which are relevant to 
the Hclfrich-Hurault instability, to zero does not lead to 
any qualitative modification of our results for the SmA- 
to-SmC transition. When Ai = A2 = A3 = 0, the model 
described in Eqs. (|5.9[) to (|5.12p is equivalent to that stud- 
ied in Ref. [26j when polarization is ignored. When c a is 
integrated out of /, the result is identical to the elastic 
energy density ffQ of Sec. |IV] with C5 renormalized to 
C5M — C5 — X 2 /(2r), with D m , m — 1,2,3, replaced by 
A m A|/r, and with E replaced by [v/A + r/2)\\/r A . 

We can now analyze the transition to the SmC phase 
in exactly the same way as we did in the strain only 
model of Sec. |lVl We complete the squares involving the 
strains and the director-strain couplings. The resulting 
elastic energy density is once more a sum of two terms, 

/ = / (1) + / (2) , 



(5.23) 



where 



= 5 C x w 2 zz + C 2 w u w zz + \ C 3 w 2 t + C 4 w 2 ab , 

(5.24) 

is quadratic in the shifted strains 



and where 

f^ = \r R c? c + \v R {<? c )\ (5.26) 

depends on c only. The coefficients o, f, and u> in 
Eqs. (|5.25p are of the same form as the coefficients o, 
t, and oj of Sec. 033 see Eqs. P~oT) and (|4~7)) . albeit with 
Di, I = 1,2,3, replaced by A;. The coefficient p is given 
by 

The renormalized elastic constants tr and vr in 
Eq. (j5T2"6l) read 



\2 



(5.28a) 



2C 5 ' 

v R = v - 2 a 2 Ci - 4 of C 2 - 2 f 2 C 3 

-2w 2 C 4 + 4p 2 C 5 . (5.28b) 

In the incompressible limit, the coefficient f vanishes 
wheras the remaining coefficients and the renormalized 
clastic constants tr and vr stay nonzero. 

The transition to the SmC phase occurs at r R = 0. 
From Table HI we have r — r t + fi(p — l) 2 /p, A4 = 
r f — fi(p 2 — l)/p, and 2C5 — r t + /i(pi) 2 /p, and we find 
that the critical value of r t at which the transition 
occurs to be zero. In other words, the coupling to the 
elastic network docs not affect the SmC transition tem- 
perature. This result is a direct consequence of the as- 
sumed semi-softness of the elastomer in the absence of 
smectic ordering. 

Next we minimize / to assess the equilibrium states. 
With our coordinate system chosen so that c aligns along 
x, we obtain readily from Eq. (|5.26[) that 



tiy = for r R > and r R < 



(5.29) 



and 



s = e 



for r R > , 

±V~ r -R/ u « for r R < • 



(5.30) 

where a is the angle that the reference-space director 
makes with the z-axis. The full reference space nematic 
director is thus 



n = (sin a, 0, cos a). 



(5.31) 



Wzz 


= u zz 


"^ 2 , 


(5.25a) 


y° xx = 1 (f 4 


■ LO — 


o)S 2 


Wu 


— Uii ~ 


-fc 2 


(5.25b) 


U° = i (f ~ 

U yy 2 \l 


- U> — 


0) s 2 


Wab 


= u ab 


- Q (c a C b - ^Sabcl) , 


(5.25c) 


u zz = d S 2 > 






Waz 


= u az 


- PCa , 


(5.25d) 


u xz u zx = 


pS , 





Note that this corresponds to a counterclockwise rotation 
through a about the y-axis of the original director n = 
(0,0,1) in the SmA phase. The director (|5.31[) is also 
the target space director under a symmetric deformation 
tensor A as shown in Fig. [5jb) . 

The components of the equilibrium strain tensor then 
follow from Eqs. (|5.26|) as 

(5.32a) 
(5.32b) 
(5.32c) 
(5.32d) 
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and zero for the remaining components. Thus, to leading 
order in the order parameter S, the equilibrium strain 
tensor has exactly the same form as the one predicted by 
the strain-only theory of Sec. HVl The only differences 
reside in the specifics of the fore-factors of the S- and 
S^-terms, which are qualitatively unimportant. 

Once again, we have to choose our coordinate system 
in target space. As in Sec. |IV] we choose this system so 
that the transition form SmA to SmC amounts to the 
simple shear shown in Fig. [IJc) with tan 4> = A° z /A zz , 
With this choice, 



and A° zx = 



A° = 



A 

XX 





A 

XZ 





A 

yy 











A 



where 



(5.33) 

(5.34a) 
(5.34b) 

(5.34c) 
(5.34d) 



Knowing c° and A° we can discuss what happens in 
the SmC phase to the layer normal, the director, and the 
uniaxial anisotropy axis. Under the simple shear (|5.34[) . 



A 

XX 




- (r + uo 


- a) S* , 




A 

yy 


= Vh 


- (r — u> 






A 

XZ 




2pS 








h(f + Q 


-a)S 2 ' 




A 

ZZ 




- 2a 2 S 2 - 


Ap 2 S 2 {\ 
1 + (f + U) 


-S 2 ) 
-a) S 2 



(A^)" 1 = (A",) -1 **, and hence 
N= (0,0,1). 



(5.35) 



Thus, as expected, the shear deformation induced by the 
transition to the Sm-C phase slides the smectic layers 
parallel to each other. In this geometry, it does not rotate 
the layer normal. Since N is parallel to the z-axis under 
simple shear, the angle between the layer normal and the 
nematic director n is the angle that the director makes 
with the z axis under simple shear. This angle is simply 
O = a+f3, where (3 is the angle through which the sample 
has to be rotated to bring the symmetric-shear configura- 
tion to the simple-shear configuration. Under symmetric 
shear, the symmetric deformation tensor is given by 



A, 



!J 



1/2 = (S + 2U) 1 ' 2 



(5.36) 



In order to calculate f3, we need the symmetric equilib- 
rium deformation tensor A° , given by Eq. (j5 . 36[) with g 

replaced by g° = A 0T A°. In terms of the components of 
P = tan" 1 ([(1 + 2 U °) 1/2 W[(1 + 2m°) 1/2 ] 2Z ) 



2pS: 



p-1 



p 



1 



s, 



(5.37) 



where we replaced r t by r£ — to obtain the final result. 
Note that u xz and (3 are positive as depicted in Fig. 



Tedious but straightforward algebra verifies that the 
simple-shear deformation tensor A , whose components 

are given by Eq. (j5~34)l . satisfies A° = £ y (/3)(g°) 1/2 , 
where 



cos (3 sin (3 

1 
— sin (3 cos (3 



(5.38) 



is the matrix for a counter-clockwise rotation about the 
y-axis, which is into the paper in Fig. O through (3. The 
angle between n and N, which is equivalent to the angle 
between n and the z-axis, is 



e = 



2p 



P 



1 



S. 



(5.39) 



The uniaxial anisotropy vector e becomes e = 
(sin/3, 0,cos/3). 

Note finally that the angle 7 in Fig.[5]is 7 = tan -1 [(l + 
2u°) 1 / 2 ] xz /[{l + 2u°) 1 / 2 ) zz w u xz . Thus, to lowest order 
7 = j3. They differ, however, at higher order in S. The 
tilt angle <f> depicted in Fig. [Tfc) is given in terms of 
the angles defined in Fig. [5]by0 = /3 + 7«2/3«2(p — 
1)/ (p+l)S. Thus, the spontaneous mechanical tilt of the 
sample, as described by <p, and the tilt of the mesogens, 
as described by 0, are not equal. 



D. Elasticity of the smectic-C phase 

To study the elastic properties of the SmC phase we 
expand the clastic energy density / about the C2h equi- 
librium state. Expansion of to harmonic order re- 
sults in 



Sf {1) = \Cx (Sw zz ) 2 + C 2 8w u 5w zz + I C 3 (Sw u ) 2 



+ C 4 {SWabf + C 5 {SWazf 

with the composite strains 



6w zz 


= Su zz 


— 2a S Sc x 


SWti - 


= Su t i - 


2f S Sc x , 


SWxx 


= ~5Wyy = \ (5U : 


5w xy 


= 5W yx 


fi^ocy 


8w xz 


= Sw zx 


= Su xz - 


Sw yz 


= SW Z y 


= SUy Z ~ , 



xx fiUyy) UjS Sc x 
U) S Sc y 

p (1 - § s 2 ) 6c x , 
p (1 - i s 2 ) Sc y . 



(5.40) 



(5.41a) 
(5.41b) 
(5.41c) 
(5.41d) 
(5.41e) 
(5.41f) 



The expansion of is particularly simple. It leads to 



Sf 



(2) _ 



v R S 2 {8c x f 



(5.42) 



A glance at Eqs. (|5T4"0|l . (|5T4T|) and ([532D shows that 
the 2 components of the c-director, 5c x and 5c y , play 
qualitatively different roles. Whereas 8c y appears only 
in the composite strains (|5.4ip . the component Sc x also 
appears in Eq. (|5.42[) . In the spirit of Landau theory 



17 



of phase transitions, the term vrS 2 (Sc x ) 2 makes 5c x a 
massive variable. 5c y , on the other hand, is massless. 
Since 8c x is massive, the softness of the SmC phase that 
we expect from what we have learned in Sec. IIVI cannot 
come from the relaxation of Sc x . Rather it has to result 
from the relaxation of Sc y . Anticipating this relaxation 
Sc y , we rearrange Sf^ so that Sc y appears only in one 
place. Then we combine the two contributions Sf^ and 
Sf^ and integrate out the massive variable Sc x . Some 
details on these steps are outlined in Appendix [Cl 

Our final step in deriving the elastic energy den- 
sity is to change from the strain variable 8u to u' = 
(A 0T ) _1 5u (A ) -1 with the equilibrium deformation ten- 
sor as given in Eqs. (|5.34[) . This takes us to 



/SmC - /c°fc + A 



6c„ 



A 



o (2ALC 4 n + A° xz C 5 E)u' xy + A°,C 5 S< ; 



A 



(5.43) 



where U = -QS,E = -p(l - S 2 /2), and A = 2C 4 n 2 + 
C5S 2 . /^? f * is exactly of the same form as the result 
stated in Eq. (|4.17[) . The only differences lie in the 
specifics of the elastic constants. Our final formulas for 
the elastic constants, which are rather lengthy, are col- 
lected in App. [B] 

Equation (|5.43[) shows clearly that Sc y can relax locally 

to 



SCy = 



-A 



(2A° XX C 4 H + A xz C 5 E) u' xy + A° ZZ C 5 E u' y: 



(5.44) 



which eliminates the dependence of the elastic energy 
density on the linear combination of strains appearing 
on the right hand side of Eq. (|5.44p . In other words, 
the relaxation of 5c y produces an elastic energy density 
identical to that of our strain-only model presented in 
Sec. HV1 



r rsoft 

/SmC — JC 2 h ' 



(5.45) 



up to the aforementioned differences in the specific details 
of the elastic constants. These details do not affect the 
elasticity qualitatively. As in Sec. IIVI the limit S — > re- 
produces the uniaxial elastic energy density of Eq. (|2.15p 
and the incompressible limit leaves the form of fsmc un ~ 
changed. Most importantly, our model with strain and 
director predicts the same softness of SmC elastomers as 
our strain-only model of Sec. IIVI 

In our analysis we have completely neglected the Frank 
energy, i.e., the effects of a non spatially homogeneous 
director. However, it is legitimate to ask if it might affect 
the softness of the material because, a priory, it is not 
impossible that the Frank energy could lead to a mass for 
Sc y . We address this question in appendix iDl where we 
find that dc y remains massless even if the Frank energy 
is included. 



E. Rotational invariance and soft extensional 
strains 

In this section we will discuss the softness of SmC elas- 
tomers from the viewpoint of rotational invariance in the 
xy-plane of the reference space. The results presented 
here depend only on symmetries and not on the detailed 
form of any free energy. They thus apply quantitatively 
even when strains are large. Moreover, we will inquire 
whether SmC elastomers are, like nematics and biaxial 
smectics, soft with respect to certain extensional strains. 
Here, we will use a somewhat different starting point than 
in Sec. IIII CI in that we first consider soft deformations 
P, 0, l5q rather than soft strains. Firstly, this is inter- 
esting in its own right. Secondly, this will set the stage 
for a comparison of our theory to the work of AW [13] on 
the softness of SmC elastomers. The results presented 
here depend only on symmetry and the existence of a 
broken-symmetry state with the symmetry of the SmC 
phase. They are not restricted to the Landau expansion 
of the free-energy we used in preceding sections or to 
small strains. 

Let us first determine the general form of soft defor- 
mations. The equilibrium or "ground state" deformation 
tensor A maps points in the reference space to points 

in the target space via R(x) = A°x. Rotational invari- 
ance about the z axis in the reference space ensures that 
R(0^^x) = A O" x describes a state with equal en- 
ergy, i.e., an alternative ground state. In other words, a 
deformation described by 



(5.46) 



has the same energy as one described by A . Any de- 
formation A relative to the original reference system can 
be expressed in terms of a deformation A' relative to the 
reference system obtained from the original reference sys- 
tem via A° through the relation 



A = A'A°. 



(5.47) 



Thus choosing A = A , we find that the deformation 



^ = A g B 1 z (A°r 1 



(5.48) 



with O „ the counterclockwise rotation matrix as given 

=-K 1 z ° 

in Eq. (|3.26[) , describes a zero-energy deformation of the 
reference state represented by A . Further rotations of A' 
in the target space, of course, do not change the energy, 
and the most general soft deformation tensor is 



A 



(5.49) 



where T is an arbitrary target-space rotation matrix. 

The reasoning just presented applies to any elastomer 
with rotational invariance about the z-axis in reference 
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space. Now let us turn to the specifics for SmC's. In- 
serting the equilibrium deformation tensor as calculated 
in Sec. lIV Al or Sec. IV CI we obtain 

(cos?? r 1 / 2 sin?? s[l — cosi?] \ 
— rj_ sini? cos?? sr^ 1 ^ 2 sin?? I > (5.50) 
1 J 

for SmC elastomers, where s = A xz /A^ z and, as be- 
fore, r± = (A° x /Ay ) 2 . Of particular interest to our dis- 
cussion of response to imposed strain, which we present 



shortly, will be soft strains with a vanishing xy compo- 
nent. To construct such a soft deformation tensor, we 
rotate through an angle u> about the z axis, 

I' = £ T ,>)A'. (5.51) 

Then, the condition A' xy = is satisfied when the target- 
space and reference space rotation angles are related via 
tan lu = r 1 / 2 tan ?? in which case 



A' 



1 -1/2 
2 r ± 



1) sin 20 





cos 2 ?? + r_i_ sin 2 ?? 




l -1/2 r 



— s[l — cos??] 
- (r± - 1 ) sin 2?? + 2r ± sin ??] 



(5.52) 



where = [1 + (r± - 1) sin 2 ??]- 1 / 2 . When ?? = tt/2, 
then 



£=&(4 = n/2) = I o rf a ) , (5-53) 

1 



-1/2 n -1/2 

r, (J —sr. 



which corresponds to an overall deformation 



A = A A = 



A 

yy 










A 

XX 


A 








A 



(5.54) 



i.e., to a shear deformation of the original Srm4 in the 
yz- rather than the xz-plane. Figure [6] shows the effect 
of deformations A for a series of values of ?? between 
and 7r/2. 

Knowing the deformation tensor A , we know how the 
shape of a sample changes in a soft deformation. An 
equally interesting question is, however, how the orienta- 
tion of the mesogens changes under such a deformation. 
To address this question, we recall from Sec. IV Al that 
rotations in the reference space do not affect the target 
space vectors such as the director n. Thus, n does not 
change in response to rotation described by 0_r j2 (??). Of 
course, under target space rotations, n does change ac- 
cording to n' i — OxAjTij. Therefore, the director, as given 
in Eq. (|5.39[) . changes in the process of the soft deforma- 
tion A to 



n= (cosu sina^, sinw sina^, cosa^)- 



(5.55) 



At this point we pause briefly to compare our find- 
ings to AW. Our soft deformation tensor (|5.51[) is, up 
to differences in notation, identical to the soft deforma- 
tion tensor found by AW, see the last equation in the 
appendix Ref. [3^|. The same holds true for the change 
in the director associated with these soft deformations, 



Eq. (|5.55[) . However, whereas the AW derivation empha- 
sizes geometric constraints, ours emphasizes that softness 
arises from invariances with respect to reference space ro- 
tations and the independent nature of reference and tar- 
get space rotations. It should be noted that unlike the 
AW derivation, ours does not impose incompressibility; 
rather the incompressibility condition of the soft defor- 
mation arises naturally from the from of Eq. (|5.48[) . 

To discuss the implication of the rotational invariance 
on the Lagrange elastic energy, we now we switch from 
deformations to strains. In terms of the soft deformation, 
the general form of the soft strain tensor is given by 



ATA' - 1 = i [(A') t A' - £ (5.56) 



independent of target-space rotations. Inserting 
Eq. (|5.50p or (|5.52|) . it is straightforward to check 
that Eq. (|5.56p and Eq. (|3.27ap are equivalent. Using 
Eq. (|530")) we find 



< K (tf) = -^-^[l-cos2??], 
u'Jd) = "-^—L sin 2??. 



(5.57a) 
(5.57b) 



u' xz (d) = — - [1 - (r_L - 1) cos??] sin 2 - , (5.57c) 

T _]_ L 



KM 



4r± 
s 

n 2 



(5.57d) 

r± — (r± — 1) cos??] sin?? , (5.57e) 



u' zz (d) = -— [1 + rx- (r_L - 1) cos??] sin 2 - (5.57f) 

Tj_ Z 



for the specifics of the soft strain. Equation (|5.57a[) im- 
plies that the soft strain has nonzero components for in- 
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FIG. 6: Effect of deformations A as given in Eq. (|5,52|l for a series of values of $ between (a) and (f) 7r/2 with $ increased 
in steps of 7r/10. In the process a parallelepiped-shaped sample with initial shear in the xz-plane is transformed into a 
parallelepiped with shear in the j/z-plane that appears as if the original parallelepiped had been rotated by n/2 about the 



finitesimal 1?, viz. 



u' xy {d) = u' yx (d) 
KM = <yW 



S 



(5.58a) 
(5.58b) 



To ensure that these infinitesimal strains do not cost elas- 
tic energy the following combination of elastic constants 
has to vanish: 



C X yXy ( r ± I) 2 + 2 C X yy Z S(r_|_ 1) + Cy Z y Z S — 



This equation is fulfilled if 



C X yxy — C COS 9 , 



Cxij 



C cos 6 sin 6 . 



Cyzyz — C sin 



with the angle 6 given by 



tan 



l-r± 



(5.59) 



(5.60a) 
(5.60b) 
(5.60c) 



(5.61) 



Note by comparing Eqs. (|5.60[) and (|4.17[) that the anal- 
ysis of the rotational invariance presented here gives ex- 
actly the same relations between the elastic constants as 
our analyses of the SmA-to-SmC phase transition pre- 
sented in Sec. |TV] and Sec. IVCltolVDl 

Modifying our arguments slightly, we can also under- 
stand from them the vanishing of the elastic constants 
C x ' y 'y>z' and Cy'z'y'z' hi the elastic energy density (|4.23p . 
Rotating the soft strain tensor with the rotation ma- 
trix (|4.22p with the rotation angle 9 given by Eq. (|5.61|) 



leads to a soft strain tensor that has in the limit of small 
i? only 



J-y'z 



(0) 



iy'z' 



(r± - l) 2 



<d (5.62) 



For this strain to cost no en- 
must be zero as they are in 



as nonzero components, 
ergy C^y > y > z > and Cy> z i % 
Eq. lEP!]) . 

Next we turn to the question whether extensional 
strains can be soft in SmC elastomers. As we did for bi- 
axial smectics we consider extensional strains along the 
y-axis, u' yy > 0, as a specific example. Again we assume, 
for the sake of the argument, positive anisotropy in the 
xy-plane, r± > 1. Equations (|5.57a|) imply that u' yy is 
converted into a zero-energy rotation through an angle $ 
as given in Eq. (|3.30[) . if the remaining components relax 
to 



-r_,_ U y y , 



u'(r± - 1 - 2< ) 



■ry 



2r_L 



2 m! 



mi 



l x 2u 'vv 



- 1 



(5.63a) 
(5.63b) 

(5.63c) 



2u 'vv 



r± (r± - 1) 



rx-J(rx-l)(r±-l-2u') 



yy 



2u 'yy 
rx-1 



(5.63d) 
(5.63e) 
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When u' yy is increased from zero to u' yy = u yy = (r± — 
l)/2, ■& grows from zero to n/2 and the state of the elas- 
tomer, originally described by the equilibrium strain ten- 
sor u° is changed without costing elastic energy to 



1 



-s/rj_ 




(5.64) 



which is, of course, the strain tensor associated with the 
deformation tensor of Eq. (|5.53|) . 

In this process, the shape of the sample changes as de- 
picted in Fig. [5] As already discussed, the configuration 
at •& = 7r/2 describes a sample in which the SmA phase 
was sheared in the yz- rather than the xz plane. Thus 
further increase in u' yy beyond u yy is equivalent to in- 
creasing u' xx beyond zero in the original sample sheared 
in the xz-plane. Thus, the second Piola-Kirchhoff stress 
is 



yy 




if U yy K Uyy, 

u yy u yy 1 



(5.65) 



where Y x is the Youngs modulus for stretching along x. 
Equation (|5.65[) implies that the stress usually measured 
in experiments, i.e. the engineering stress, is given at 
leading order by 



yy 




Kyy) 



if a: 

if A: 



T J ^II (5.66) 

yy 



> 



Therefore, when plotted as a function of the deforma- 
tion A' yy , the engineering stress a yy s for a SmC elastomer 
looks qualitatively the same as the corresponding curve 
for a nematic or a biaxial smectic elastomer, cf. Fig. [31 



VI. CONCLUDING REMARKS 



an elongation in the y direction (the direction perpendic- 
ular to the order in the smectic plane) as a specific ex- 
ample. However, we could impose an extensional strain 
in any direction and we would find softness, provided 
the c-director has freedom to rotate into that direction. 
This excludes only stretches in directions lying in the 
plane spanned by the equilibrium c-director and the ini- 
tial uniaxial direction (i.e., in our convention stretches in 
directions lying in the xz-plane). Of course, the width 
of the soft plateau in the stress-strain curve, cf. Fig. [3J 
depends on how much the c-director can rotate until it 
has reached the direction of a stretch. Therefore, the soft 
plateau will be most pronounced for stretches perpendic- 
ular to the plane spanned by the equilibrium c-director 
and the initial uniaxial direction (our y-direction) . 

The softness of SmC elastomers is more intricate than 
that of biaxial smectics. At first sight it seems as if it 
takes a very specific combination of shears to achieve a 
soft response. However, with the coordinate system cho- 
sen appropriately, it turns out that SmC elastomers are 
soft with respect to certain conventional shear strains 
(with our conventions shears in the y'z'-plane). Even 
more important, as far as possible experimental realiza- 
tions of softness SmC elastomers are concerned, is that 
these materials are also soft under extensional strains. 
What we have said above for biaxial smectics also ap- 
plies here: the experimental boundary conditions must 
be right and the direction of the imposed stretch must 
be so that the c-director can rotate. 

As pointed out in the introduction, very recently con- 
siderable experimental progress was made by Hiraoka et 
al. , who synthesized a monodomain sample of a SmC 
elastomer forming spontaneously from a SmA phase upon 
cooling. This is exactly the type of elastomer for which 
our SmC theory was made. Therefore, it seems well 
founded to hope that our predictions for SmC elastomers 
can be tested experimentally in the near future. 



In summary, we have presented models for transitions 
from uniaxial Srm4 elastomers to biaxial and SmC elas- 
tomers: Landau-like phenomenological models as func- 
tions of the Cauchy-Saint-Laurent strain tensor for both 
the transitions as well as a detailed model for the transi- 
tion from the Srm4 to the smectic-C phase. The detailed 
model includes contributions from the elastic network, 
smectic layer compression, and coupling of the Frank di- 
rector to the smectic layer normal, and allowed for esti- 
mating the magnitudes of its phenomenological coupling 
constants. We employed the three models to investigate 
the nature of the soft elasticity, required by symmetry, 
of monodomain samples of the biaxial and SmC phases. 

We learned that biaxial smectic elastomers are soft 
with respect to shears in the smectic plane. In addition 
to that we saw that extensional strains can be converted 
by the material into zero-energy rotations, provided the 
experimental boundary conditions are not too restrictive 
and allow the remaining strain degrees of freedom to re- 
lax. We illustrated this softness by explicitly considering 
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APPENDIX A: EFFECTS OF HIGHER ORDERS 
IN THE STRAINS 

A priori, a formulation of stretching energy densities in 
terms of the variables r\ and r\ z provides a more adequate 
framework for discussing the incompressible limit than a 
formulation in terms of the respective linearized expres- 
sions uu and u zz . Using ua and u zz , on the other hand, 
makes our models more tractable and perhaps also more 
intuitive because we stay in close contact to the standard 
formulation of elasticity as presented in textbooks. The 
purpose of this appendix is to discuss what changes occur 
in our theories if we use 77 and r\ z instead of uu and u zz . 
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Biaxial elastomers 



for the soft biaxial state, where 



Now we use Eq. (|2. 12|) as the starting point for setting 
up our model for soft biaxial elastomers. When C4 can 
become negative, higher order terms featuring u ab , rj and 
r) z must be added to Eq. (|2. 12|) so that the model elastic 
energy density becomes 



f (2) = f 
J uni J u 



+ A 1 t 1z u 2 ab + A 2 r]u 2 ab + B{u 2 ab ) 2 



(Al) 



Proceeding in close analogy to the steps described fol- 
lowing Eq. (|3.ip , we find that the equilibrium values of 77 
and r\ z are 



2 ' 



(A2a) 
(A2b) 



with a and /3 as given in Eq. Eq. (|3.5[) . The equilibrium 
values vP az and v? ab remain unchanged. To learn more 
about the equilibrium state, our next task is to determine 
the equilibrium strain tensor u . Using our knowledge 
about u° az and it is clear that u° is of the form 

?(t + S) 

?(t-S) I , (A3) 











where we used the abbreviation t 



t and u„, are 



unknown thus far and we need to determine them as 
functions of the order parameter S. Replacing r\ z and 
u zz in Eq. (|2.11bj) by rp z and u° zz , using Eq. (|A2a|) for 77° 
and solving for vP zz we find 



(1 



■ s 2 r-i 



(A4) 



To leading order in S this result reduces to vP zz = \aS 2 
which coincides with the result of our linearized theory 
presented in Sec. It remains to calculate t. Since u° 
is diagonal, Eq. (|2.11a[) leads simply to 

rj° = [(l + t + S)(l+t-S)(l + 2u° zz )} 1/2 - 



Using Eqs. 
t: 



t = 



(A5) 

2]) for rf and u° zz we obtain by solving for 



/3 5 2 )2 



(1 + f S 2 ) 5 



S 2 



1/2 



1 . 



(A6) 



To leading order in S this expression reduces to t = \{j3— 
a + 1)S 2 . Comparing this to the result vP cc = |(/3 — a)S 2 
of our linearized theory, we see that the linearized theory 
lacks the contribution ^S 2 to v? cc . Nevertheless, up to 
this detail, the form of the equilibrium state predicted by 
the linearized and the non-linearized theory is the same. 
In the end, the non-linearized theory leads to the elastic 
energy density 

/ =\C l9zz 2 (5u zz ) 2 + C 2 g-^u8u zz + W 3 (u) 2 



+ Cr (8u az ) 2 + B 2 S 2 (6u xx - Suy y y 
+ {AiSg-^ 2 Su„ + A 2 Su)(5u xx - 



SUyy), (A7) 



= [dct(l + 2 1q )] 1 / 2 U+ 2 io )- 1 5 i ] 4 



Vb , sy_ 

, ^ (A8) 

where Vq is the volume of the reference uniaxial state, 
Vb is the volume of the biaxial state and SV — V — Vb- 
Thus, volume changes in the biaxial phase are suppressed 
by the term proportional to C3 in Eq. (|A7[) . In the in- 
compressible limit C3 — > 00, the nonlinear theory indeed 
produces fixed SV — and not Sua = as the linearized 
theory does. Our findings about the softness of the biax- 
ial state, however, remain practically the same. 



2. SmC-elastomers 

When C5 can become negative in response to SmC 
ordering of the mesogens, higher order terms featuring 
u az , T] and rj z must be added to Eq. (|2 . 1 2[) to ensure me- 
chanical stability. Then the model elastic energy density 
becomes 



f (2) = f 
J uni J 1 



uni -r Dx r\ z u az 
+ E{u 2 az f. 



D 2 Tj U az + D 3 UabUazU bz 

(A9) 



Essentially repeating the algebra described following 
Eq. P~Tj) . we find 



r lz > = a(u G ] z ) 2 =aS 2 , 
°-T(u° az ) 2 = TS 2 , 



V 



with a and r as given in Eq. (|4.6 



(AlOa) 
(AlOb) 

and u„, remain 



the same as in Sec. IIV A[ Consequently, the equilibrium 
strain tensor is of the form 



(All) 




where t — 



For determining the equilibrium state, it 



remains to compute t and u° zz as functions of S. From 
Eqs. (|2.11b|) and (|A10a|) it follows readily that 



1 r 

2 



fl + ffS 2 



(A12) 



which is to leading order in S identical to the result 
vP zz — tr S 2 of our linearized theories for SmC elastomers. 
Equations (|2.11bj) and (IA11|) lead us to 



V °=[(l + t) 2 (l + 2u a zz )-4(l + t)S 2 ] 1/2 -l. (A13) 



Taking into account the equilibrium values (|A10a[) and 
solving for t we find 



t = 



2S 2 



(1 + aS 2 ) 2 



AS 4 



(1 + tS 2 ) 2 

{1 + aS 2 ) 2 {l + aS 2 Y 



1/2 



- 1 . 
(A14) 
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To leading order in S, t — (r — a + 2) S 2 which has to 
be compared to the vP cc = (r — a) S 2 stemming from 
our linearized theories for SmC elastomers. Note that 
the linearized theories miss the contribution 2S 2 to v? cc . 
However, as was the case for the biaxial soft state, this 
discrepancy does not affect the nature of the softness. 



APPENDIX B: ELASTIC CONSTANTS 

This appendix collects our results for the elastic con- 
stants of biaxial and SmC elastomers. 



1. Elastic constants of biaxial elastomers 

The elastic constants of soft biaxial elastomers as de- 
fined in the elastic energy density (|3.23|) are given to 
order 0(S 3 ) by 

C zzzz = d + 2C 2 + C 3 + 2a[d + 2C 2 + C 3 ] S 2 , (Bl) 

C xzxz = 2[C 5 + C 5 S+ \{a + /3)C 5 S 2 ] , (B2) 

C yzyz = 2[C 5 ~C 5 S+\{a + 0)C B S 2 } , (B3) 
C zzxx = C 2 + C 3 + [C 2 +C 3 +Ax+ A 2 ] S 

+ [\{ a + P){C 2 +C 3 ) + A 1+ A 2 ]S 2 , (B4) 
C zzyy = C 2 + C 3 - [C 2 + C 3 + A 1+ An] S 

+ [^(a + (3)(C 2 +C 3 ) + A 1 +A 2 }S 2 , (B5) 
C XX xx — Cs + 2[C 3 + A 2 ] S 

+ [(1 - a + (3)C 3 + 4A 2 + 2B] S 2 , (B6) 

Cyyyy = C3 — 2[C 3 + A 2 ] S 

+ [(1 - a + P)C 3 + AA 2 + 2B] S 2 , (B7) 
C xxyy = C 3 -[(l+a- 0)C 3 + 2B] S 2 . (B8) 



For the remaining elastic constants we find to order 
0(S 3 ) 



C zzzz =Ci + 2C 2 + C 3 + 4(a - 2)(d + 2C 2 + C 3 ) S 2 , 

(Bll) 

C xzxz = 8[2(d + 2C 2 + C 3 ) + 2{D 1 + D 2 ) + E]S 2 , 

(B12) 

C zzxx = C 2 + C 3 + [{a + t + u - 4)(C 2 + C 3 ) 

+ A(C 1 + 2C 2 + C 3 ) + A{D 1 +D 2 )]S 2 , (B13) 

C zzvv = C 2 + C 3 + (* + t-lo- 4)(C 2 + C 3 ) S 2 , 

(B14) 

C xxxx = C 3 + C 4 + 2[4C* 2 + (4 + t + u - <j)C 3 

+ (r + u ~ a)C A + 2(2D 2 + D 3 )] S 2 , (B15) 

Cyyyy = C 3 + C 4 + 2(t - U ~ (j)(C 3 + C 4 ) , (B16) 

C xxyy = C 3 - C 4 + 2[2(C 2 + C 3 ) -(a- r)(C 3 - C 4 ) 

+ 2D 2 - D 3 ] S 2 , (B17) 
C xxxz = [4(C 2 + C 3 ) + 2D 2 +D 3 ]S, (B18) 
C vyxz = [4(C 2 + C 3 ) + 2D 2 - D 3 ] S , (B19) 
C zzxz = 2[2(d + 2C 2 + C 3 ) +D 1 +D 2 ]S. (B20) 



b. Elastic constants as obtained from the model with strain 
and director of Sec. [V| 



For 9 and C we find 



2. Elastic constants of SmC elastomers 

Here we list our results for the elastic constants of soft 
SmC elastomers as defined in Eq. (|4.17|) . 

a. Elastic constants as obtained from the strain-only model 
of Sec. \TV\ 

Our results for the angle 9 and the elastic constant C 
read / _»o tt \ 

6 = tan_1 An - Z \n tt = S + °(^ 3 ) , (B21) 

- tan Ul^Ai^sy ( } ' ( } c-i C4Ce ( A0 „) 2 ( A ^ S - A ^ n ) 2 

(Ag,) 2 (AL-Ag^^) 2 A , cos l d 

C = ^ ^0 -iq | ^ + [2^-^ + 2^^ ^ 

= 4C 4 + 4C 4 [2(t - a) + (w - 4)w] S" 2 + 0(S* 4 ) . ^ 2 C s 

(BIO) +0(5 4 ). (B22) 
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The remaining elastic constants are given to order 0(S 3 ) 
by 



APPENDIX C: STEPS LEADING TO EQ. (f5^3)) 



c z 



C x 



C1+2C2 + C3 + 4 (a-2p 2 ) 



x (d + 2C 2 + C 3 )-2 



\ 2 

A 4 



cv, 



S 2 , (B23) 



2p 2 (d+2C 2 + C 3 ) 



3 C 5 + A 4 [A 4 -2(2A 2 -A 3 )] 



A 2 



5'- 



(B24) 



CzZ£C£C — C2 + C3 



4p 2 Ci + (4p 2 + a + f + Q)C 2 



+ (a + f + u)C 3 - 4 



2A 2 : + A 1 (2A 2 + A 3 ) 



C 5 



PA3A4 - A 2 (A 3 + 2pA 4 ) 

+ 4 — L-5 



A 2 



(B25) 



C Z zyy — C 2 + C3 + 



(a + f - uj - 4p 2 )(C 2 + C 3 ) 



i (A 1 + A 2 )(2A 2 -A 3 ) Ca 



Ai 



C X xxx — C3 + C4 + 2 



(B26) 



4p 2 (C 2 + C 3 ) 



+ (f + Q-a)(C 3 + C 4 ) 

(2A 2 + A 3 ) 2 + 8pA 4 (Ai +A 2 ) 



X 2 
A 4 



S' 2 



C 3 + C 4 + 2 
(2A 2 - A 3 ) 2 



\2 

A 4 



(r-a)-a)(C 3 +C 4 ) 
S 2 , 



(B27) 



(B28) 



Cxxyy — C*3 C*4 + 2 



2p 2 (C 2 + C 3 )- (<7-f)(C 3 -C 4 ) 



(A 4 + A 2 )(2A 2 - A 3 ) + pA 4 (2A 2 + A 3 ) 
^ — 05 



\2 
A 4 



Cxxxz 4 

Cyyxz 2 

Czzx.z 2 



p(c 2 + c 3 ) - — 



Ai + A 2 



s, 



2p(C 2 + C 3 )- 2 -^±^C 5 
M 



2p{d + 2C 2 + C 3 ) - 2X \ h C 5 



(B29) 
(B30) 

S, (B31) 

,5'. 



(B32) 



In this appendix we outline some of the algebraic steps 
leading from Eqs. ([5H)f and (j5T42j) to Eq. lfQ5| . As 
discussed in the text following Eq. (|5.26p , the role of 8c v 
is special in that it is the local relaxation of this quantity 
that makes SmC elastomers soft. To see this, we recast 
Eq. (|5.40p . where Sc y appears in two different terms, by 
using 



2C 4 (Sw X y) 2 + C 5 (Suiyz) 2 



■■ A 
2 



SCy 



A 



2C 4 Id Su xy + C 5 S Suyz 



[TISUyz - ESUxyY 



(CI) 



with IT, H, and A as defined below Eq. (|5.43[) . The va- 
lidity of Eq. (jCip can be checked by straightforward but 
slightly tedious algebra. 



The second major step in going from Eqs. (I5.40|l and 
(|5.42p . after these are combined, to Eq. (|5.43|) is to inte- 
grate out the massive variable 8c x , i.e., to replace 8c x by 
its minimum value. This minimization is straightforward 
but a little inconvenient because of the number of terms 
that are involved. Combining Eq. (|C1[) and the outcome 
of this minimization we arrive at the intermediate result 
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Sf = A 
1 

+ 2 
1 

+ 2 



2C4 II (5m™ + C5 ^ Su yz 



S£y + 



C 1 +2C 2 + C 3 --(X 1 +\ 2 ) 2 S 2 



C3 + Ci 



9. 



A 2 + 



S 1 



{8u 



f 2 * 5 [II 5u yz - S 5u xy ] 2 
(6u zz f + 

1 



C 3 + C 4 + — ( A 



2 

S 2 



(SUyyf 



c 2 


+ c 3 - 




hA 2 ) 


(a 2 + 


A 3 ^ 
2y 


S 2 




c 3 


-c 4 - 




hA 2 ) 


(a 2 - 


A 3 ^ 
2y 


S 2 





C 2 + C 3 --(A 1 +A 2 ) (^A 2 - 
-^(Ai + A 2 )A 4 s(\-^S 2 



A3 
2 



8U ZZ 5Uyy 



\aS\1 



8Uyy8UxZ 



1 



A, 



A, 



A 4 5( 



S 2 

8u zz 8u xz , (C2) 



where we have used the shorthand £1 = gS* 2 + C5p 2 (l + 5S 2 /2). Switching from the strain variable Su to it' = 
(A 0T ) _1 %(A°) _1 then takes us to Eq. (j5^gj) . 

I 



APPENDIX D: EFFECTS OF THE FRANK 
ENERGY ON THE SOFTNESS 

In Sec. |V]we have neglected the Frank energy. In this 
appendix we check if our conclusions about the softness of 
SmC elastomers are affected if we take the Frank energy, 
or rather the corresponding density 

jWnk = \K X [V • n] 2 + \K 2 [n • (V x n)] 2 

+ \K Z [n x (V x n)] 2 , (Dl) 

into account. To this end we have to check if Eq. (|Dip 
can lead to a mass for Sc y . 

Each of the three terms in Eq. (|D1|) features a deriva- 
tive. Since the director n lives in the target space, 
these derivatives are derivatives with respect to the tar- 
get space coordinate R, i.e., 

Vkni = d£ = ^^di- (D2) 

In Sec. [V] we modeled SmC elastomers via converting the 
director to a reference space vector by using the polar 
decomposition theorem. Doing the same conversion here 
we get with the help of Eq. (|5.3p 

^(A-^j^+^SO' <D3) 



The second term in the braces contains a derivative of 
hj and hence it is clear that it cannot lead to a massive 
term. This leaves us with the first term in the braces as 
a potential source of a mass. Using Eq. ()5.2|1 we get 



= ^VA ij +-, (D4) 



where rjAij — \{djUi — diUj) is the antisymmetric part 
of the displacement gradient tensor 77^ = djUi. Thus, 
the first term in the braces couples fij to derivatives of 
rjij and higher order terms, and as such this term cannot 
make Sc y massive in a perturbative (diagrammatic) ex- 
pansion. Furthermore, the terms in this expansion will be 
subdominant compared to the u yz c y at small wavenum- 
ber q. Rotational invariance in the reference space should 
dictate that Frank energy terms do not ever generate a 
mass for c y , a full demonstration of this fact to all orders 
in perturbation theory is beyond the scope of this paper. 
We conclude that our findings about the soft elasticity 
of SmC elastomers remain unchanged if we include the 
Frank energy in our model. 
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